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ABSTRACT 


Pressure Distributions, Added-Mass, and Damping Coefficients 
for Cylinders Oscillating in a Free Surface 


A linearized theory is developed for the pressure distribu- 
tions, added-mass and damping coefficients for horizontal cylinders 
oscillating vertically with small amplitude while semi-immersed in 
the free surface of a fluid of uniform depth. The results of sam- 
ple calculations made by digital computer are presented. The cal- 
culations of the pressure distribution and the total vertical 
force for the case of a circular cylinder in fluid of infinite 


depth are compared with values measured in an experiment. 


The problem is formulated as a linearized boundary-value 
problem in the theory of infinitesimal surface waves. A velocity 
potential is synthesized from appropriate functions. The boundary 
conditions on the free surface and on the bottom and that the sur- 
face takes the form of outgoing waves far from the body are satis- 
fied exactly. The boundary condition on the surface of the cylin-e 
der is satisfied at the rest position of the cylinder by an expan- 
sion in non-orthogonal functions. Convergence in the general case 
is assumed. 


The expansion coefficients for the case of a circular cylin- 
der in fluid of infinite depth were calculated by two different nu- 
merical. procedures. One of these methods was chosen for further 
development to provide sample calculations for elliptic and more 
general cylinders of approximately ship-like cross-section. 


The pressure fluctuation at several locations on the surface 
of a circular cylinder and the total vertical force required to sus~ 
tain vertical oscillations were measured. Two-dimensional conditions 
were simulated. The measured values and values based on the sample 
calculations show similar distinctive behavior. 


The results of a general study of instrumentation for the 
measurement of pressure fluctuations in similar experiments are pre~ 
sented. 
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The study of tha nydrodynamics of cylinders oscillating ina 
free surface is of interest for several reasons and has therefore 
received attention for some years. Ursell [3949] formulated the 
problem of the horizontal circular cylinder, semi-submerged and 
undergoing forcec periodic heaving oscillations, as a troundary-value 
problem in the theory of infinitesimal gravity waves. He completed 
the rigorous solution in principle and presented the results of a 
numerical calculation that leads to added=mass and damping coeffi- 
Clients. One might say the modern history of the problem begins with 


this solution by Ursell. 


Grim [1953] follows a somewhat different procedure in synthe- 
Sizing a potential soluticn but with a restricted number of coeffi- 
cients so that the results are approximate. He compares his circular- 
cylinder results with Ursell’s calculation. Grim extended his method 
to present approximate results for certain elliptic cylinders and for 
certain cylinders of somewhat ship-like cross-section called "Lewis 


Forms" in the field of naval architecture. 


fasai 11959] returned to the method of Ursell, using in prin- 
ciple an unlimited number of expansion coeffisients, and, by a mapping 
procedure, extended the results of Ursell to include elliptic and 
Lewis-form cyiinders. He presents the results of numerical calcula- 
tions for the added-mass coefficient and waveheight ratio (equivalent 
to damping factor) for certain elliptic and Lewis-form cylinders. 
Tasai compares his results for an elliptic cylinder (Chorizontal-to- 


vertical-axes ratic of 1,5} to the results of Grim. Grim’s approximate 


* References identified in the text by author and date are collected 
at the end. 





results are markedly iow for the added-mass coefficient. There is 


closer agreement for waveheight ratio. 


The work of Ursell, Grim, and Tasai cited above is restricted 
to water of infinite depth. Yu 12960] added the condition of uniform 
finite-depth to the probiem as formniated by Ursell. The potential 
functions used by Yu are those appropriate to infinitesimal gravity- 
waves in water of finite depth. Yu completes the solution in princi- 
ple following the methed of Urseil for a circular cylinder. He pre- 
sents the results of a numerical calculation for the added-mass coeffi- 
cient and for the waveheight ratio. His result recovers the solution 
of Ursell as the water deptn approaches infinity; practically speak- 
ing, the numerical results are essentially the same when the water 


depth is equal to ten radius-lengths. 


Yu also reported an experiment with a semi-submerged cylinder 
in simulated two-dimensional couditions.e He measured the waveheight 
ratio and shows substantial agreement with his calculated results. 
The few other reported experiments simulating two-dimensional condi- 
tions are not readily compared. Holstein £1936] used a prismatic 
cylinder of somewhat shipelike cross-section, nearly rectangular. 
Schuler [1936] used a similar form but in free oscillation, a very 
different hydrodynamic problem. Dimpker [1934] used several pris- 
matic forms, including a circular wylinder, but again reports mea- 


surements for free oscillaticne, 


There exists work on the similar problem in three dimensions, 


on approximate methods for solving the three-dimensional problem by 


drawing on twoedimensional resvits, end experimental work with models 





not simulating two-dimersion comuditiens. dicivding the third dimen— 
Sion is a sufficient departure from the objective of this work that 


pursuit of that worthy problem is not undertaken here. 


The present paper considers the case cz the more general cylin- 
ders in a fluid that may be of uniform finite depth. The problem is 
formulated as a boundary-value problem in the theory of infinitesimal 
gravity waves. Within this framework the problem is linear and the 
solution is obtained by superposition cf appropriate velocity poten- 
tials. Following Urseil, the soluticn is obtained in principle employ- 
ing an unlimited number of expansion coefficients. For the special 
case of the circular cylinder, one recovers the solution of Yu. For 
the special case of infinite water depth, a solution is obtained which 
includes that of Tasai and is also explicit for more general cylinders 
than he treated (for example, the three-parameter cross-section family 
of Landweber and Macagno L1959]). Finally, for the special case of 
infinite depth and a circurar cylinder, the parent solution of Ursell 


is recoverec, 


The hydredynamic pressure distribution for the cases of the 
previcus soluticns can pe derived from the work of the authors cited, 
However, none state this result er present sample calculations for 
their cases. The pressure cistribution is considered here in order 
to provide possible insignt inte similar problems of naval architectur- 
al interest and to provide calculations for comparison with experiment- 


al measurements. 


The results of numerical calculations for the case of a circular 


cylinder in fluid of infizite depth are presented. The numerical 





procedure differs from that of Ursell for the same case. The force 
and waveheight ratio results are in agrsement. The voressure distribu- 


tion was not previously reported. 


{f} 


In order to compare the theoretical results with experimental 
values, measurements were made on a circular cylinder in simulated 
two-dimensicnal conditions. Total vertical fcrce and water pressure 
fluctuation on the model surface were measured, The measured values 
compare favorably with the predicted values. Over a small part of the 
frequency range used, a distinct departure from two-dimensional condi- 
tions was observed and was due to the particular experimental arrange- 


ment. In this small range, measured values consistently deviate from 


predicted two~-dimeltisional values. 


The results of numerical calculations for elliptic cylinders 
and certain cylinders of ship-like cross-section in fluid of infinite 
depth are presented. Certain added=-mass ana waveheight ratios of this 
group can be compared with those of Tasai, who used a different numeri- 
cal procedure. The pressure distributions and the results for certain 


cross-sections have not been previously reported. 


Numerical calculations following this procedure for circular, 
elliptic and more general cylinders in fluid of finite depth are to be 
done at the Computer Center, University cf California. The calculations 


reported here were done at this facility. 





II. Formulation of the Problem 


rea Description. 


We consider cylinders of a certain class to be defined more 
precisely in the next section. The cylinder is immersed in a pre- 
viously-undisturbed fluid which may be of uniform finite or infinite 
depth. The axis of the cylinder is in the plane of the still free 
surface. Consider now that the cylinder is forced in vertical sim- 
ple harmonic motion and that steady-state conditions are attained, 
We assume two-dimensional conditions, and that the amplitude of 
oscillation is small compared with the cylinder beam. In the steady 
state, the forced vertical motion generates a surface disturbance 
that at a distance from the cylinder takes the form of uniform pro- 
gressive gravity waves whose amplitude is proportional to the oscil- 


lation amplitude. 


We require the pressure distribution on the surface of the 
cylinder, the vertical force required to sustain the oscillations, 
and the ratio of the waveheight far from the body to the oscillation 


amplitude. 


The total vertical force in the steady-state is periodic and 
includes two components of particular interest. Because of steady- 
state periodic conditions it is possible to resolve the vertical 
force into components in phase with the motion acceleration and 
velocity. These components do not depend on the previous history 
of motion, given steady-state. I+ is possible therefore to define 
added-mass and damping coefficients as ratios of the appropriate 


force component to the mass of the displaced fluid that are functions 





of body geometry and oscillation frequency only. Expressions for 


the added-mass and damping coefficients are required. 


Ese The Cylinders. 


The cylinders admitted are those mapped by conformal trans- 


formation of circles by 


N 
Z =(2ne2) 
a G a 2 font] C 


n=0 0 
) 
-iQ 
where* a ss 
er be i pe : 
: . =i B 
%Z= xX + iy = ire 
@) . 
The elementary case ee 6) : 
provides a circular cylinder. The case B 
N=0O yields elliptic cylinders. Lewis Jy 


[1929] introduced the use of this transform 


with N=1 for a two-parameter family with ass a selected to 


y 


provide cylinders of more or less ship-like cross-section. The 
solution of Tasai 11959] employs this two-parameter family. Landweber 
and Macagno [1959] use this transform with N = 2 and select ays az 


a to provide a more general selection of ship-like cross-sections. 


2 


Prohaska 11947] treated sections derived by transformations employ- 


a and 


7 Sr aes 


a and a 


ing a 5 ae 


In this work we leave N unspecified in the solution of the 


problem, and select various a5 for sample numerical calculations. 


n+ 1 


*The n- and y-axes are chosen positive down from the free surface 
and the angular coordinates 90, B are measured from the vertical 
to facilitate reference to related work cited. 





(ae, Infinitesimal Surface Wave Theory. 


The preblem is solved within the framework of infinitesimal 
surface wave theory. We may apply this theory if we assume that 
the fluid is incompressibie, that the flow is irrotational, and 
that the motion amplitudes and velocities are all sufficiently small 
to reject all but linear terms. fn this preblem we also neglect sur- 


face tension. 


The process of linearization is specifically employed three 
times: First, in linearization cf the free-surface boundary condi- 
tion; second, in satisfying the boundary condition on the surface of 
the cylinder at the rest-position; third, in the approximation of 
the hydrodynamic pressure by using the linearized form of Bernoulli’s 
equation. These are weli-establisned elements of the general theory 


within which we work. For a more complete discussion see Wehausen 


L1960] or Stoker [1957]. 


264 The Boundary-Value Problem in Potential Theory. 


Given the assumptions ahove, w2 cam state the boundary-value 
problem in potential theory. We seek a velocity potential 9 which 
° ¢ 4 2 ° a ® 
is a solution of the Laplace equation Vy = 0 in the fluid domain 


and which satisfies the following boundary conditions, 


(i) The linearized free-surface condition: 


62 
D 
ze: 


on the surface y = G outside the body of the cylinder. Here 


K = w/Z6 





(ii) The normal velocity at the lower boundary of the fluid 


vanishes: 


on the bottom y = h. 


(iii) The radiation condition: The disturbed surface takes 
the form of regular progressive outgoing gravity waves at large 


distance from the cylinder. 


(iv) The normal velocity of the fluid at the surface of 
the cylinder is equal to the component of the cylinder’s forced 


Mewecaty in that direction. 





LII. Solution cf the Proble 


ell Descri pero. 


The prooslem has been stated in Chapter II. The solution is 
described in this section and outlined in this chapter. The com- 
plete velocity potential 9» (and the corresponding stream function 
y) describing the fluid motion will be synthesized as the sum of 
appropriate potential functions. Each component function individu- 
ally satisfies the Laplace equation and the boundary conditions on 
the free surface and on the bottom. The asymptotic behavior of the 
component functions for large distance from the cylinder correctly 
contributes to a potential appropriate to outgoing surface waves, 
and thus the radiation condition is satisfied. In order to satisfy 
the last boundary condition, concerning the normal velocity of the 
fluid at the cylinder boundary, an expansion in noneorthogonal 
functions is required. It is shown that an unlimited number of 


coefficients ca 


Ny 


5 ih principle, be determined. Convergence of the 
expansion, proved for the circular cylinder by Ursell (1949, 1953!, 
is assumed for the general case. Thus the complete pocvential and 
stream function are known anc the problem is solved. Having the 
potential, we determine the pressure distribution and, by integra- 
tion, the vertical force required to sustain the oscillation. The 


vertical force is resolved into components in phase with the ac- 


09 
(J. 
4 
1) 
9 


celeraticn and velocicy to dded-mass and damping coefficients. 
The ratio cf waveheiga*t +o oscillation amplitude is available from 


the asymptotic behavior (progressive surface waves) of the potential 


and from the camping ccefficient. 








10 


Dae Notation of the Geometry. 


The z= or physicai-plane geometry is conformally mapped from 


the reference C-plane by the transform 


Zo =(2nt1) 
a =¢ ye Banei S 
n=O 
where a and Bonny are real so that the ze-plane coordinates mapped 
19 


Prom ¢, = ipe™* are 


a 
Z=x+ iy = ire F 9 





; s a aes ‘ 
x= al psin o + 2(-1) P| sin(2n+1)0], 
n=0 p 
N a 
y = al pcos @ + Ga) aS oo cos(2n+1)0| ° 


n=0 pao 


Special notation will be convenient for coordinate values on the 


cylinder surface- mapped from the reference circle — = pee 


riWOGy J) = ake. eo X(p os @) 
x(body, 9@= n/2) = aG, ca= X(P os t/2) 
y(body) = ak, i Ypos Q) 


The product aG is the half-beam, b, of the cylinder at 


the free surface. The fluid depth is uniform at y =h. 


SAPs) The Component Potential Functions. 


In the following sections we shall develop the complete veloci- 


ty potential » and stream function y in the form 





W 
oO £ en Fon P05 
sin «wt + ece wh+ 7 P>,( Os wt + 4 Sin wt | + 
“f wal ¥ 
i = "Om irs 
WwW 
ct ‘D %s 
a om “2m 
+ q. : iE ’ 
eee ron a ON ae a ut = WwW cos wt ° 
2m "2m 
Components are named as follows: 
ve source potential; 
2. standing-wave potential; 
Pon? dom expansion coefficients; 
Po multipole potentials 
w : | , — 
P50 associated standing-wave potential. 
° Pe 
The combined components sin wt + -_ cos wt and the combined 
s S 


Ol COS a. Lorn sin wt 


components a : Pe, ee display asymptotic behavior 
¥;,, Sin wt Y5, cos wt 
2m em 
Pom 
appropriate to cutgoing surface waves. The multipole an turn 
Vom 
consists of two parts, ¢ath with several terms. 
Be 1. The Multipole Potential Pom ° 
The muivipole potential is written in two parts 
a rey i) f = % its ¢ Y = 2 z 
Fer Contd. + Lomi Xo¥ ot) m 10 9/9°°0° 


t' (r ,B) satisfies the free-surface condition and 
the bottom ccndition fer infinite depth but does not satisfy the 


bottom cenditiocn fer finite dapth. The part $55,,(Xo¥ oh) is added 





2 


to make the total satisfy the finite-depth boundary condition, follow- 
ing Thorne [1S52)6" The vfunction Pom (Xs ¥ oh) independently satisfies 
the free-surface condition and for h +a, 059 ¥s¥ 9h) >O. 


The index 2m will be understood when it is dropped from the 


symbols Pom and Pom ° 


' 
(a) Part Pon Fo 8B): 





mir. a) = aa 2_2mp mets 2 ql emt ent 2) 55 (2m-2n+2) 6 
2n+1 Om+2n+2 7 
n=O r 
: K a~"cos(2m-1) 8 K ye -1)2 een act+ ent? 45 (2m+2n+1)B 
2me=1 pomeaL é 2m+on+L pout 2n+ 1 ae: 


To show that y'(r,B) satisfies the free-surface condition, 


it is convenient to express that boundary condition in another form: 





0” \ 
om 7 j d” " dg"! 2 a a 


The part Ky' + d9'/dy , independently zero, is written 


Ky! — = $e Onn GO mam 2.5 


When the substitution is made and evaluated at B = n/2, the remain- 


ing terms are 


ee mt + Bienen eee oe + y(a1)%a nee [cos sn+ 
=O 2n+1 


& BE Ee ey ee 


We have introduced 2s = 2m+e2n+2. Now, using the identities 

m m+1 : 
cos mn = (-1) and sin(2m-1)n/2 = (-1) , it is clear this sum 
is identically zero; therefore, g'(r,B) satisfies the free-surface 


condition, 





ake, 


The function 05 n6F9 8) satisfies the bottom condition only 


in the case of infinite depth. The total a 5 + P50 will be 


m 


adjusted to satisfy the bottom condition at finite depth in (c) 


below. 
(b) Part Poni Xo¥ sh) s 
oD 
@ (x,y,h) =/ lc, (k)sinh ky + c,(k)cosh k(h-y)]cos kx dk, 
O 


where c,(k), c, (x) are found below and assure convergence. 


Carrying out the operation (K + @/dy)p leads to the 
definition 


k 
OVS) 2 a) Eas st Se era 


to assure that o'"(x,y,h) satisfies the free-surface condition 


ay. ¥Y = Oe 


rc) The Finite-Depth Boundary Condition for the Multipole. 


The coefficient c_(k) in the integral term o'(x,y,h) 


gy, 


is determined by using the bottom condition, 


oo am i 
ey Gy dy 


To find the first term, use the expansion (for y > 0) 


cD 


i 
= ~~ ke 
n (n=1) 3 7 


Dee Oe aie dk . 


[Whittaker and Watson, A Course of Modern Analysis, Cambridge 1927. 


para. 12.2 Ex.2]. With this expansion, one may write for g', when 


yoo, 
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’ ° a : 
© = / Lote coskxy ak 
Q 
eMmeo2 Ou, 4 N ee aoe 4c 
ae a” (K+k) a a (K+k) 
where Lk) ee ee eee 2. \eok) a 
\eme1l)} Rate ( eT : 


The bottom condition is satisfied when 





: em ag ok 
aed Bue ~ + oj(k) k cosh kh } coskx dk 
Oy u & 
eel O 
vanishes. The definition c,(k) = tes Jeosh kh assures this. 


With c(k) determined, @5(Xs¥ ph) is written as a Cauchy 


Principal Value integral (). 





a 8) 

= -kh K sinh ky - k cosh kv 

P5(X,y¥ sh) a Le pa Bn Sy = cOSn TY coskx dk 

: ar Pa ¢.¢5'3) aaa -2 2m ee ee _ 

ae )) | K cosh kh = k cosh Kh ~~~ : 
Somes 
1 ye) a é ies 9) a i K sinh ky - k cosh k ae 
n=O ioral: (25—1)) K cosh kh - k cosh kh °° 


Bic The Gombinec Multipole and Associated Standing-Wave Potentials. 


The radiation condition [2.%%iii)] requires that the asymptotic 
behavior of the potential for large x and y= 0 be appropriate to 
of the maltipcole 


outwardly progressing surface waves. The part 


‘ 4 * ° 
he part © of the multipole 





ct 


tends to zero for |x; ->@. However, 


does not tend to zers for |x| Pp, y = 0. The associated standing- 


a W ; ae 
wave potential 95, is shown below to combine with the multipole to 


form a total which satisfies the radiation condition. 


@ » om’ = o,h) @p"(x,0,h) for 


fy 
© 
ry 
ys 
ef 


| x | =>@. It is shown in Appendix Ae-l that the following asymptotic 


expansion hoids: 





aby 


cosh K (hey) 


tt e 
~5(x,y¥ 5h) ~~ = Eon “cosh Kh Sin K |x| as | x | >O, 
and 
1 ; 
G5(X,0,h) aE sink. |x| as |x| > OD, 
where K em em 
Ree CKnenke hh) eee | 0 : 
2 Oe 0 eK oh + sinh ek h (2m-1)! 
N a5 1 K ae 2s 
na (2s-1)! 


and K \(K, bh) is defined as the positive real root of 


K=K tanh Kh. 
O O 


For finite h >O the asymptotic behavior of ''(x,o0,h) 
as | x | >om with an appropriate harmonic function of time, is 
identified as typical of a standing wave in a fluid of depth h. 


We add a standing-wave potential consisting of the function 
Ww 
Po, = So, Cosh K (h-y) cos K x malye. > iene 


multiplied by a harmonic function of time. The coefficients Con 


are chosen in the following so that the combined potential repre- 


sents progressive waves and the radiation condition is satisfied. 


To determine ¢ consider 


em’ 
WwW " WwW , 
PonC0S wt + 95 Sin sere ss 05_(%s0sh)cos wt + Q5,,(%,0,h)sin wtas |x| >@. 


The right-hand side is asymptotically equal to 


= sin K |x|coswt+ ¢ cosh kK hcosK x sin wt as |x| >. 
O em O O 


om 
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By defining 


the asymptotics benayior of the 


ac 


combined functions is 


. 4 he roy \ 
- Bin (Kix) = We) , 


4 


which displays the required behavior. 


We also form the sum with conjugate harmonic time functions, 


P sin wt - 
‘ 2m rene iy OS 
which for | x | 2?>am, y = O, tends asymptotically to 


- EA, cos (Ki |x| = wt) . 


A linear combination of multipole potentials, adjusted as 


above to be progressive-wave producing, is 


oD oD 
(oar ae a 
2 Pen (95,008 wt + ee sin wt) + 2. dan Pam Sth wt Pom cos wt) o 


The coefficients p, (Ka) and Qo, Ka) are to be determined, 


and 
are assumed to behave so that convergence is assured, 
The asymptotic form of the free surface y=0O for large 
me >O 1s 


poem | 
It 
09 JY 


GD 
ee aa: COS K. Ke wt ) = = > oe EB sin (K) X=wt) o 
B mel 8 m=) 


In complex notation, if we reserve the operator i 


for 
geometric variables ana introduce the operator 


of time 


j for functions 
%, thie sum is 





ib 


aD 
f a ° f r , W o | wt 
SF age Pom * 9 2m? Pom * J Pom)® : 


| 
M== 4, 


We will later make use of the corresponding stream functions 


Ww + . 
You and Yom? which may be found from 


( = { 7 4 
P5462) ‘One om 
of 1)™,o0,-<™ “ eee ( 1)8,728 1)” ik aon =(2m—iy) . 
— = ‘Onell es Pmel zs 
n=O 
- ; (-1)"a. _(21)8 —tK 228, -Ces-1) _ 

n=0 en+1 2sel 
a ogee sinkz+ k cos kz a 


= 3 K cosh kh ~- k sinh kh 


Ww W Ww 
and f, fb 2) a Aas ab Yon ° 


Cc 


ao cosh K (heiaz). 


With these definitions, the complete wave-producing multipole is 


OG 
~ W i, = ie 
Bem (ose) 3 a5) )2o4Ga) + 3) fy (ze. 


Soy, The Source Potential. 


The potential function that represents a source at the origin 


and that satisfies the free-surface boundary condition is 


oO ky _ ao _-mx 
€ ESOS Ke dk = ea err a le (m cos my - K sinm dm . 
K-k a 2 
s) m+ K 


It is easily seen from tne form on the right that this satisfies 





18 


the free-surface condition. The bottom condition, however, is not 
satisfied in the case cf fiuid of finite depth. To satisfy the 
finite-depth bottom condition, we add a term in the manner of 


Thorne [1953], 


@ 
if Le, (Kk) sinh ky + c,(k)cosh k(he-y)] coskx dk. 
O 

The coefficients c,(k), c5(k) are determined below and assure 


convergence, 


This integral term satisfies the freeesurface condition 
given that 


k 
Bei) = ioc, (k) oes eda eh 


The coefficient c, (k) is determined so that the combined 


term 


. K cosh k(h-y) 


K cosh khek sinh kh 


-ky aD 
E cos kx ) ' 
Kok dk + ¢ c, (k) (sinh ky = 


0 .@) 


cos kx dk 


satisfies the finite-depth boundary condition. This boundary 


condition, @p/dy = 0, y = h, leads to 


e L 
ey) = Gare coer” 


With this coefficient, the source potential is 


= cosh k( hey) cos kx 


vas a K cosh ki = k sinh kh 


dk . 


5) 


When h?q@, this simplifies to the source potential 


given first. 
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3.34% The Combined Scurce and Starding-Wave Potentials. 


The radiation condition is satisfied by adding a standing- 
wave potential to the sources potential, with the necessary harmonic 


functions of time, to form the term 
cos We + © sin wt 
ape w a Wo o 
The standing-wave potential function is 


~, = a, cosh K Chey) cosK x y 


with ay defined as given below, 


Appendix A=e2 shows that the asymptotic behavior of the source- 


function ». for Ix] ~#@, y=O is 


B_ sink |x| 
8 O 


where en{ cosh K h)* 
B (Kb! eK bh + sinh 2K h ° 


Now a. is defined as 


S 
o cosh Kh 


a 


so that the asymptotic form, as |x| PFO , 


(» coswt+ y_ sinws) wa coshkK bcos Kk, |x| cos wt + 
C Ss hoe Oo O 
y= 


+ B sink |x| sin wt 
= 0 


may be written 


by 
3 
Oo 
{fq 
ms 
an 

Ss 
o 


td) 
“; 
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The source potential and its associated standing-wave potential 


OD & 
cosh k(h-y) cos kx ; S 7 
i ao eee eee Kh cosh K (h-y) cos K x cos wt 


therefore satisfy the free-surface, finite-depth, and radiation bound- 


eaey CONdItLiONS. 


In the special case h ™qm@m this may be written 


~Ky ar aed cos my - K sinmy , ~-Ky 
ne sinkxe /{ 5 5 an | Sin wt + neéE 
O m + K 


cos Kx cos wt , 


which is the form used by Ursell [1949] and Tasai [1959]. 


mhe form Of the frée surface for large “xX > O for ether 
case is 


WE 


7 = boon 
& 





© sin (Kx ~ wt) . 


SY 


For the special case h>@, Es el Se 


The stream function corresponding to this potential may be 


found from 


Re, (£4(z) + j£,(2))e 9% 
where 
pez =P se TY 
en cosh Koh 
= 3K h + sinh 2k n °OS K Ch+iz) 
f(z) one ee 


a cosh K,(h+iz) 


= € X cosh kh- & sinhkh °* ° 





rial 


34 The Potential Solution. 


The potential solution written as a function of the complex 
variable z= x +iy is 
aay = =" Re Ct (2) + jf (2) + Z (p +Jqa5,/(f, cz tits me 2) de ge 
TW i ee oe S art om om 
The coefficients Po Kas qs,,( Ka) are determined later by use of 
the normal-velocity boundary condition. The right-hand side has 


been normalized by the coefficient J(K oh) in order to display the 


surface waveheight hy at large distance from the cylinder. 


The function J(K oh) may be determined as follows, for exam- 
ple, from the asymptotic development of the potential for large 


map O, y= O. 


gh 


. WwW 
OE ol 5S ( >) ty to fy ce s a 
ee, o,h) ~ =A [B.co K xe) + z Pon® oy Sin (XK x-wt ) 


es dono a Cas (K x=ut)| 
me] ‘ 


The corresponding asymptotic surface waveheight as x *a@ is 


oils) & =()% 
= (E. - z fame e ay sin CK. x=wt) + 


+ Dis te cos (K xewt) | ° 
This may be written 


hse CK x -]-wt oF 
- i (Ko (OG cr v) 





when we cefine 


wD OD 
c ne oe | 2 st 2 
& { ; = = | 4 ant a hy ; * ) f \ 
Me: : Nt eae 
oD uD 
, = 
tan t= 2 Pobag/ By ~ 2 donP an) * 
mei m=1 


For 8h ?@, we recall that EL, 20 , E Sr, so that J #1. 


Ali quantities in the potential and stream function 


W 
ty 0) = 
A gh, aS. ee eee oe Pom . om 
=. =a Sin ws + sos we + ei COS Wu + Sin wt ‘f 
TG : 
Y Ys Vs m= Vom Yom 
go on ; - ea : ty 
a a Dou, ( | Sin ad wW cos W J 9 
at Vom Vor = 


except the expansion coefficients for a given cylinder, water- 
depth, and heaving oscillation have been determined. When the 
expansion coefficients are determined in the next section it will 
be convenient to have the following notation, The substript a 


designates evaluation on the cylinder boundary. 





gn. 
Ld ig anne nal : “ 4p 
oJ = ‘M sinat+ N cos wt) 
: ku!) 
boay 
CD GM 
(kez 5 = 5 (5 fe 5 O.. fa 
Baa Vea ke 2 Pom “ema a 2 Goma 
m=]. a | 
© x 2 
N a = =- dp = 
(Kaye) OCR f dom Pome 2 Bem oma: 
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gh 


W 
y J See (D sinwt+ C coswt) 
body 


aD >) 
W 
D(Ka,@) = Ysa * 2 Pom Yoma * 2 dom Yoma 


m=1 
aD cD " 
C(Ka,@) = oa * a Pom Yo ayo tom YQma ° 
m=1 ma m=1 
In addition we designate the evaluation of y] for O=B8 = nf2 by 
body 
gh 
yJ eg (B sinwt+ A coswt) 
body fe 


oD aD 
W 
B(Ka) = D(Ka,n/2) = Ysa n/2* 22m Yoma n/2* 282m Yoma 1/2 


oD cD 


W 
A(Ka) = C(Ka,n/2) = UCR aye 2 Pom Yoma 1/2 7 22m Yoma n/2 . 


The apparent redundancy of notation has the purpose of facilitating 


reference to related solutions previously mentioned. 


BoD Evaluating the Expansion Coefficients. 


The vertical velocity component that is perpendicular to 
the cylinder surface must equal the normal velocity of the fluid 
at that point. This boundary condition is to be satisfied at the 
rest positicn of the cylinder by suitable choice of expansion co- 
(Ka). 


efficients (Ka), 


Pom tom 
The vertical motion is y = h) cos (wt+ ¢). After using 


the substitutions 
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Ox _ oy _ , 
oe a 8Y5 t 
a 98 bLody 
the boundary condition is a 
n 
_ 4 _ 2 ox) y V 
Qs dt 0s ° 


body 


We evaluate this expression at the special location 8B = 6 = n/2, 


using the notation of the preceding section: 


gh 


yJ = a (A cos wt+ B sinwt) = - ai ag F 
body,n/2 " 
Then we may use 
gh 
a = a (A coswt+ B sinwt) 


to evaluate the same expression in general locations on the cylinder: 


y J =a (C cos wt + Dsinwt)= - Lx) 
body e body 
gh 


> Ge (A coswt+ B sinwt)aXx. 


Equating the coefficients of sinwt and of coswt, we obtain 


0) 


be 
Lad 

i) 

Oo 


C(Ka,@) - A(Ka) 


a 


X(0) 


~ 
© 

il 

O 


D(Ka,9) - B( Ka) : 


a 





ea, 


Substitution of the definitions of A, B, C, and D and introduc- 


tion, for convenience, of the notation 


QP 


£5 KAO) Yoman/2 ~ Y2ma ° ee nee 


Qy>s 


W W 
Bo \ Kas 0) Voman/2 ~ Yoma ? Maly eyyeees 


leads to the two equations 
aD aD x 
et 7. ea era 2. dom52m ~ Yea ~ G Yean/2 ? 
m2 ice 


cD aD 
Geen Se a Aon ton =) Vea a Jes 
m=1 m=1 nee 


These two equations display the roles of Po, (ka) and d5,, (Ka) 
as expansion coefficients. Ina given example the coefficients 
can in principle be evaluated. In practical calculations some 
Suiticient finite number of coefficients can be found by the ap- 
proximation methods employed for expansions in non-orthogonal 


functions, as in the sample calculations given later. 


In the special case of fluid of infinite depth, (Ka,0) 


bom 


is identically zero. 


B10 The Ratio of Waveheight to Oscillation Amplitude. 


We return to the expression for the normal velocity condi- 


tion evaluated at B = @ = nf/2 given above. 


gh 


W : 
a, (A coswt+ B sinwt) 


dy 
- st (aG) 


wh \b San (Wt wie one 
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The trigonometric expansion for sin(wt+ ¢e€) leads to 


ane 6: 
wb oo : 
A=--m=—-—- ndJsine ., 
g h 
W 
co. oh 
@W bo 
Bzsmas=—--—-~ nt JCOSE . 
gs h 


= 


These give, respectively, by division and addition after squaring, 


Ew Kb nt J 
4G in + B 


The first result defines the previously-unspecified phase angle 
€é of the motion. The second gives the ratio of the height at 
infinity of the generated progressive waves to the oscillation 
amplitude. 


507 Pressure on the Cylinder Surface. 


Gauge pressure at a point on the cylinder surface with 


reference to the constant pressure P, on the free surface is 


(P = p.) = pay - p $e 


from the linearized form of Bernoulli’s equation. The first term 
on the right contributes a hydrostatic pressure due to the loca- 


tion of the point at rest, 
yo BL , 


and a fluctuation of the hydrostatic pressure due to vertical motion 


pgh, cos (Wwte+ €) o 





a 


e7 


The hycrodynamic pressure is 


Fy 
{ 


Q 


pe wy 


= [= M coswt+ N sinwt] 





where we have used the notation for the potential o» introduced 
in Section 3.4% and M‘(Ka,@), N(Ka,@) are to be evaiuated at the 


Maine on the cylinder of interest. 


The acceleration cf the given motion 


y =h. cos fwt+ ce) 


h 
— (B coswt- A sinwt) 


: Ae b Bo 


2 
| Crh. . 
it =) ; 
jy O = ==eeee=== (A sinwte B coswt) . 


ta? 4 ar 


yp is resolved into components in phase 





ay 


Tne hydrodynamic pressvur 
feeietcne acceleration ani wlth the vertical velocity: 


i* 


Pes MB» NA 


pKa 6.t) = a 5 (A sinwt = B ces wt) + 
6 ode Fs os R 


“” 
hemgieMihe SO NB ® ol to 
<> =e aS . afi COS 60% < B S231 “re ) i] 
= L- + 2 


After substituting for the waveheight ns this may be written 





>= et MB + NA y i pb MA = NB ae 
——————— a, = > ry 
| A” 4 Be ie + ce 


Srein units of the hydrostatic cluctuation peh, . 
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p aaa MB + WA yay Ky & 
ogh ar Aveo 
fers A Ba RA A + 5B 


4 
hota 
= 
Av ive 
on 


where M(Ka,@), N(Ka,@) ars ewaluated for the location of interest. 


The total charge in fluid pressure must include the fluctuation in 


us ae) 


hydrostatic pressure due co vertical motion. This may be added in 
normallized form y/h, to the above expression to give the total. 


Examples for specific cylinders are given later. 


Bice The Vertical Force Required to Sustain Oscillations. 


One part of the verticai force is that opposing the result- 
ant of the hydrodynamic pressure. This is called the hydrodynamic 
force. The hydrodynamic force per unit cylinder length is found 
by integrating the vertical component of the hydrodynamic pressure 


acting on an element of area. 


rey) 


0° 


ta 


amo = pl Ka pect) eo 


S o 


After substituting cosa= 6x/fds and also : 


- N a (2ntijas 4 
BE. ds = dx =a P cos OQ + of = oe COS (2061) | de 
= aW(6)de 
in the form 
WCQ) 


aGxX«= b —=G ce 


where >» = a@ is the haifebeam at the sree surface, the force 


becomes 





eg 


ne ray 
r= 2b ff pi(Ka,9,t) a 
0 


a@ . 


Recalling the notation for the pressure introduced in the last sec- 


tion, one may express this in the form 





epgh. b 
1 oe 2 (N  sinws= M cos wt) 
KX oo e) 2 
wrerse 
yes 
Rd pe ON eD 
Mt Ka) = f Mi Ka,@Q) ta oe 5 
o 
{2 
ey . 
_ WOE) 
NA‘ Ka) a ah N“ Ka,@) ~a ea) ° 
cy = 


The hydrodynamic force is resolved following the procedure 


Semene preceding section, 





apeh 2 M je x 
Ke a 4 A 
F(Ka,t) = aay SS = (A sinwt = B cos wt) 
ee Com 
AU + B 


eth bo M AN B 
, W we) wD 





ee 8 COS ot B Bae 
a a* " ae 


substituting for the waveheight D9 che finds 


“Be NA 
‘4d 


cee me gis . ee Bee Mom niet 
FE = i is oe aS 5 dn =p Ky 2 eel US? ° 


A » Bo A + -&B 


It i8 convenient So norma lize the hydrodyeami: force in units 
of the change of hydrostatic force (naval architectural "tons- 


4 
pereinch immersion"), 








a a" i ea 2. oy ee 
eegbh, ~ Ko 3 2 ee ) + Kb —=5 Bt mae ee 
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The total force fluctuation is found bv adding the normalized 

buoyancy change v/h and of course account must be taken of 
9 J rot 

the inertia of the mass of the body and net buoyancy at the rest 


position. 


B.9 The Added-Mass and Damping Coefficients. 


Since the fluid motion, and thus the pressure on the cylin- 
der surface and the vertical force, is periodic, it is meaningful 
to define addedemass and camping coefficients and they will not 


be functions of time, i.e. of the past history of the motion. 


3e91 The Added Mass and Added-Mass Coefficient. 


The added mass is the ratio of the hydrodynamic force in 
phase with the acceleration to the acceleration 


MB+NA 
fe) O 


2 2 ss 


ye 
m A’ +B 


Tl 
ma 5 b p 
The asymptotic value of the added mass for large Kb in the spe- 


cial case of infinitely deep fluid is given in Section 4.4. 


It is convenient in this problem to define an added-mass 
coefficient as the ratio of the added mass m to the mass of the 


fluid displaced by a circular cylinder of equal beam. 





M B + NA 
m 4 o O 
qZ.c2 ft 2 2 ‘ 
5 p A’ + B 


To normalize the added mass with units of the actual displaced 


fluid mass, one may use tne volume of the cylinder per unit length, 





pee 


) 
, ne on ® \ Kae 
(2n+1, Bong? 


(2nsl)2.. =| ? 
= we) 


so that the ratio of the added mass to the mass of the displaced 








fad is 
2, MB+NA 
ae 6 a : 
p \ Tt Ac Be 


The coefficient defined in this manner is inconvenient for the 
case of a fiat plate oriented perpendicular te the motion. In 


this case V= 03; Rmowever, the added mass is still defined and 


finite. 


5 09a The Radiated Energy and the Damping Coefficient. 


An energy balance must exict between the average rate of 


‘a 


£3 


work done in one cycle of vertical motion amd the average energy 


fiiux radiated by the outwardegoing surface waves. The average 


= 


rate of work done over cone period T is 
4 a tr 
eng yf 
ey) 


Where c¢ is the factcr of proportionality for the comporent of 


hyarodynamic force in phase with the velocity, 
Ea -eCoy at) * 


The average energy flux per unit width of periodic cutgoing waves 


on two sides is 


Ce 


ee a 
ae oe 





De 


Therefore, 


2 





h e 2 
c= (-*) AB = mw pp Bow. 
O GQ) Av+B 


When this is compared with the results in Section 3.8, 
2a 
MA - NO Baad /2 o 


This is useful for checking numerical calculations and simplifies 
the expression for the component of vertical force in phase with 


the motion velocity, 





eee 
F, = 2p" Sed SS wy 
A” +B 
The ratio 
= z ond’ wy 
aoe Oe ee 
5° Pp At 5 


leads to the convenient use of 


onde 


ac + Bo 


as a damping-force coefficient. If the added-mass coefficient is re- 
ferred to the actual mass of fluid displaced by the volume V, the 


corresponding definition of the damping-force coefficient is 


ro 
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IV. The Special Case of Infinite Fluid-Depth 


The case of a fluid of infinite depth deserves special men- 
tion. An important reason is that the expansion to satisfy the 
normai-velocity boundary condition on the cylinder can be done with 
different functions. A second purpose is to summarize the simpli- 
fied definiticns. Finally, with respect to the added-mass coeffi- 
CGient, a correction factor accounting for free surface effects if 
readily defined. For practical calculations, this case is simpler 
and is of interest since the infinite-depth results are approached 


rapidly with increasing depth ratio. 


= Special Definition of the Multipole Potential Pom ; 


The multipole potential which may be used in case h = @ 


is designated by Pom and defined by 


aT ht 


( (2m+en+1) 


N 
2 Tie 1 : 
_ £Os 2mO | ae )e Ge 
n=O 2m+2n+1) p 


es : 
om Sa (2m-1)p°*" a 
This replaces PonT 98) and since Pom >0O as hom, there need 
be no counterpart to Papi Xo¥ ok) It is evident that Pom vanishes 
for |x| +a and thus has no role in satisfying the radiation bound- 
ary condition. It is consistent that om = QO for this case, (ft 1s 


Shown below that Doan satisfies the free-surface condition. 


whe free=ssurtace condition 


Ky + Take OF, @=+ 0/2 , | >| Boal) 


can be written 
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= 


: Res 
Kr + ee eee podlia — O 
ay - i ns ran — 
. 67/30: Ow 


+ 1/2. With 


by noting that @x/00 = C when 6 


(2n+l)a 


N 
ov n+l 2Y . 
ae = “lp sin@+ >) (#1) Saal Be ean (2n+1)9| 
n=0 
a iP i EBT eet “mT | 1 OF Ene 
n=O fp 
the free-surface condition is 
: 9=+n/2,|x| > db. 


Ka s (2n+1) Sere ees 
oe 2 a Smet |? * 397° 
n=0 p 


Simple substitution paralleling Section 3.31(a) demonstrates that 


Pon 9 like PoniF9B) s Satisfies the free-surface condition for all 
My, Ne 
The corresponding stream function is 
~ _ sin2md | Ke | S223 (2m=-1)0 ; $ 21) (éntl)a,..7 sin (2m+2n+1)0 
Yon ~ em ; (2meq1 ) ‘ (2m+2n+l ) (2m+en+1) 
p (2m=1)p n=O p 


The stream function evaluated on the cylinder is used in satisfy- 





ing the normal velocity condition. The part i j is designated 
= body 
Yona? 25 before, and at 6 = n/2, 
7 = Ka(-1)%2 1 . Sasa ipe FN 
2ma,n/2 ons Omeol 4 , (2m+2ntl)| ° 
(2m-1) p. n=0 (2m-2nt1) p, 
4.2 The Source Potential and the Complete Fotential Solution. 


The definiticn of the combined source and standing wave 





BS) 


» 


potential Co. coswt+ ©. sinwt) remains as in Sections 3.33, Bisse 


where it was written for the special case as 


K = mx Om cos m K sinm 
c YT sin kx « f anninae ser sue sin wt + 
Oo m + K 


@Ky 


+ TE cos Kx cos Wt o 


The conjugate stream function Y. sin wt + Y, cos wt is 


2 


OD oN X ’ 
-ne*¥cos Kx+ f &—(a sinmy+ K cos my) en Seer 
fe) m + K 


ne "AY sin Kx cos wt . 


For this case the potential and stream functions equivalent to 


those given in Secticn 3.4 are 


ep her, 4 ~ a) ~ 
_ WwW 2m cS} em 
a) & + Fas ) cosut + ( + 2 ee ) stn] 


C mz] er ‘ie =] - You 


The expansion coefficients Pon KB) Io, Ka) are to be evaluated 


below following the procedure of Section 3.5. 


4,3 Evaluating the Expansion Coefficients. 


The normalevelocity boundary condition on the cylinder 
surface is used to evaluate the expansion coefficients. Since 


¥Cp.) = - Cdy/dt)x( pods the two equations 


C(Ka,@) = Ae. A(Ka) = 0 , 
D(Ka,o) - S42 B(Ka) = 0 





36 


are found as before, but now with 


as = Vean/2 ‘ 2 Pen Yoman/2 ~ Gain) Na 
Wa = . 
f Ah este es ' 
OS aa “san/2 ~ 232m Yoman/2 ~ D(Kayn/2) , 
WD f= ~—] 
eka 70? — ‘oa * 2, Pen Y2ma ° 
@ —= 
D{Ka,@) = Ysa * 2. 42m Voma ° 
m=). 


With these definitions, the expansion equations are 


CO 


> Po, fo, =¥.. = =e Y 
oy 2m “2m ca G can/e2 
a dee fa =: ¥ XO) | 
<> = 9 
az om 2m ea G san/2 
where fon is now given by 

= X(O) _— - 

‘ Tread wwe o 
fami Ka@) = G  ‘2man/2 Yoma ‘K89®) ° 


This is the set of expansion functions which may be used in the 


infinite-depth case to replace the equivalent unbarred set. 


Completion of the problem and the results are identical 


in form with the corresponding steps in the general case, 


44 Asymptotic Value for the Added Mass and Definition of the 


FreeeSurface Correction. 


The added mass is the same as in Section 3.91, 


DURES 2, allropeon 
~ 2 rom 2 2 


A +B 
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The asymptctic value of the aided mass for large Kb (in 


this case of infiniteiy deep fluid 


) 


Neves? 


is 


2 eC 


2 q. ( -L)f¢ | 
= = Lp a “2 :  poealendaer 125 
— f) : =—r a, SS et ee en 
@ 5 a2 (Fo a) nad , enelje : 


which may be found in a manner following Landweber and Macagno [1957]. 
We recall that the origin of the result called mo here is in the 
work of Lewis 11929]. The ratio of the asymptotic value m, %o the 
mass of the fiuwid displaced by a circular cylinder of equal bean, 
nb“p/2, is the inertia-coefficient C of Lewis [1929] or the added- 


mass coefficient C, of Landweber and Macagno (1957, 1959]. 


The ratio of the added mass m to its asymptotic value am 
defines a free-surface correction factor that is a function of the 


non-dimensional freguency Kb only. 





The value of no for an elliptic cylinder is nop /2, from the 
definition above. Therefore, in the special case of elliptic cylin- 
ders (including the circular cylinder), the free-surface correction 
factor Ky 4s numerically aqual to the addedemass coefficient de- 
fined as the ratio of the added mass to the mass of the fluid dis- 


placed by a circular cylinder of equal beam. 


Examples of the free-surface correction factor Ky were 
Calculated for various cylinders and are presented in Part B of this 


work, 
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The approach used to define Ky could be used to define a 
correction factor for the free-surface effect on the added mass of 
cylinders in fiuid of finite depth but the asymptotic value ns 
for the case of finite depth is not yet simply expressed. 

A first-order approximation of the added mass m for ellip- 


tic cylinders valid for Kb =O is given by Ursell £1949] for the 


infinite-depth case: 


hele: 2 a 
m = = pb i= log. Kb = log. (1 + 5) + 00235] ° 


Here d is the draft, y(B=0), and »b remains the half-bean, 
x(B=n/2). It appears from various calculations that this is a 
useful approximation for other than elliptic cylinders, but this 


is not established as a theorem, 
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PART B. Numerical Analysis and Results of Sample Calculations. 


Vo The GCirctuiar Cyiinder in Infinitely Deep Fluid 


Ths principle objective of the first sample calculations was 
to provide the caicuilated pressure distribution on the surface ofa 
Beeeular cylinder in fluid of infinite depth. In addition, the 
addéed-mass and damping cosificients were calculated. These coeffi-e 
cients were calculated by Ursell and the present calculation con- 
firms his results and extends the range of the frequency variable. 


Our method of calculating the expansion coefficients Po, Ka) s d> (Ka) 


and of finding the added mass coefficients is notably different. 


The calculations were done on the IBM 704 digital computer 


at the Computer Center, University of California, Berkeley. 
p Cr, 2 9 


Bs J. Calculation of the Expansion Coefficients. 
po 1 Choice of Methcd. 


The definitions of most of the functions required in the cal- 
Culation of this case are simplified. The mapping is trivial: all 
a are zero, and pa = 1 is convenient in that the unit circle 

en+1 a) 
in the reference plane corresponds to a cylinder ef radius a. The 


angular coordinates are equal, 8 = @. The specially defined Pom 


of Chapter IV is identical with the general case oon of Chapter III. 


We will use 6 for the frequency parameter (here Ka) in all 


numerical calculaticns. 


The expansicn equations of Section 4.3 become in thie case 





t 


QD 
D P>,,(6) £6559) 


Te Yon (O99) = sinO y,,(5,2/2) 


fe) 
Ltn’ £6598) Yoq 6529) - sind Vgq 6 O0t/2) : 


where 


X(@) = sinO®, Y(@) = cose, 


Mitje) =Ge=l, 





Yon 6699) Brien (Oona sin (6 sin®) , 
-5 cos @ = a 
Ysa 6619) = € L(y + log_5 + 2. nim ¢°5 n@) sin (6 sing) - 
oO ,n 
(9+ 2) ——sinnO) cos (6sino9)], 
ey nin 


m-l 


£6658) = = sin2m9 - j Csin (2m-1)0 - (-1) 


iy sino]. 
We have substituted a form of Y¥5q6 599) from Appendix A-3.4. A 
method must be chosen to find P56 8) s a>,66) from the expansion 
equations for m=1,2,3,.00, their computation being terminated in 
practice at some finite m. The first method chosen was to generate 
m simultaneous linear equations each for Po ,(m=1,2,+6-m) and 

Qo (ML y2y 000 ym) by proceeding as follows. Each expansion equation 
was multiplied by sinkO@dQ@ and cos j@dQ@ with k=1,2,3,... and 
jeO,1,2,..- terminating at k, (j-1) so that m=k+j. The m 


equations (each for were integrated over the interval 


Pom? Ion) 
0<0<n/2. This will be called the "method of integration." The 
m equations were solved for m values Pon? Ion by the Gauss 
elimination procedure. Yu [1960] used the method of integration 
but with cosr@, r=1,2,3,.0.+,m as the only multiplier in his 


calculations. 





by 


The second method used for sample caiculations, was to gener- 
ate m simultaneous linear equations from each expansion equation 
by evaluating each at more-or-less arbitrarily chosen angles eo. rn 
the range O < oy < nf2 with 4=1,2)350009M o Urselli {1949] and 
Tasai [1960] chose to take values M=1 2 Fy 0009m, with m,< m So 
that m values cf Pon? Ion could be found from m simultaneous 
equations by the method of least squares to satisfy the truncated 
expansion equations at 9 = e. in the least-squares sense. In the 
sample calculations reported here the m equations were solved di- 


rectly by the Gauss elimination process for m expansion coeffi- 


cients. This will be called the "discrete-angle" method. 


In all cases the calculated expansion coefficients were tested 
by substituting them in the left-hand side of the basic expansion 
equations (Section 5.11) and comparing this result with the calculated 
right-hand side at various values of 6 = QO, » This calculation was 
printed from the computer with five significant figures. In the case 
of the calculations by the discrete-angle method, test values 9 = 9, 
that equal the discrete-angles e, lead to exact agreement and gave 
proof the matrix was properly inverted. For interstitial O19 the 
difference in right=- and left-hand sides is an accuracy index not 
generally zero. This index would not, in general, be zero for any 


0O< 6, < n/2 in the case of expansion coefficients calculated by the 


eo 


method of integration from integrated equations and substituted into 


the basic expansion equations. 


An example of this index for various calculations at 6 =1 


and ©, = 20°, 80° follows. In these calculations by the discrete- 


angle method the nearest ©; differed from e, by at least five degrees. 
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IRHS = LHS| 


The tabular entry is RES x 100% 

Method: Discretew-angle. Integration. 

m: 6 9 6 8 

for Po pit) s 

eo, = omg 14.52% 0.045 0.49 0.24 
80° 0.75 0.51 0.83 0.89 

for Gon 6t): 

Q, = 20° | 4.37%  Ool3 1.38 0.60 
80° 1.6 aaa dl 7Or A Oatte 


Another index used to test the calculations is reported 


Geoection 5.22. 
5 ca Wa Form of the Equations for Calculation. 


The expansion equations used in the method of integra- 


tion are as follows: 


iebg Pon: 

m n/2 sin k@ n/2 sin ko 
2. poo) f oC @) do = f Y..(6, 8) de 
m=1 a O en cos je O oa cos j9@ 

n/e sin ko 

- f yy... (6,n/2)sine ag 
ca 
O cos j9 


Ka) 0235 gas oy 
4204 ij 2 eaten, Ce) 


eae 


B | 





ue, 


for don: 
m n/2 sin k6 n/2 sin ko 
ie to: ? £ rf . ane 
2 do,8) J £5,,06,8) _ o6 = f  ¥5,(848) ae - 

m=1 Q cos je O cos jO 

T/e sin kQ 

- f v (6,n/2) sind dQ 
sa 
Oo cos j@ 


Keine, ee eee 
dO, U2 ee Ge) 
m=k + Z ° 
All integrations in these equations were reduced to the form 
of rapidly converging series or to algebraic results, before doing 


machine calculations. The forms used are given in Appendix A-3; the 


integrals on the right-hand side are not widely published. 

The expansion equations used in the discrete-angle method are 
given in Section 5.11, above. 
Doe The Principal Dependent Variables. 

5.21 The Results for AC&), Bd). 


The dependent variables A(6&) and B(6) are calculated from: 


cars 


m 
AC6) = 2 po (b)¥55,659%/2) + ¥Q_68s0/2) 
moi 
m 9 ymat 
2 6 > mane Dag OP + wsin 6 
re 
m F) 
B(6) = 2 a5,68) Yomalbsn/2) + Yon 65o%/2) 


}-' 


— ja Li. 
= 2. pe 


_—2 i 


a5,68) = [+ Si(5)] cos (6)+ Ci(6) sin(6). 





Lyd 


The substitution for Yous 50%,/2) employs the identity 


op =X 

ere *) & ? % e 
ge) ccna) cece ein) 
5 6 +k 


where Si(S) and Cit6) are the sine and cosine integrals and are 


available in tabulated or series form, 


The quantities A(5) and B(&) in the form ne + ce are use- 
ful for interpolation of the added-mass and damping coefficients, as 
pointed cut by Urseli [1957]. The phase angle of the motion 
= ang B(5)/A(&) is smoothly varying and nearly linear in 6. The 
results of this calculation and that of Ursell [1949] for A(&), BC) 


are given graphically in Figures 1 and 2. 


5.22 The Quantities M(5,@), N(6,@), and M(8)s N (5) 6 


The dependent variables M(5,@) and N(6&,@) are required for 
the calculation of the hydrcedynamic pressure distribution. They may 


be calculated from theixv definitions in Section 3.4, which in this 


case besome: 


% fry (5.0) ¢ ox 
M(6,@> Com’ O/ Pama? @) + Pua 66 Q) 


i 
iMey 


L 


5 
(cos 2m + mor °o8 (2me1)0] + 


t 
MB 


eo) 
. ‘om 
L 


i 


m 


n 
4 70 COB Orig, 2. =a sinn®) sin(6sin0)+ (y + log_6 + 


rake 
eal 
. ,. ma £08 n@) cos (5 sin@)] 9 
mo” 
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O Ursell [1949] 
+ This work 


FIGURE 1. THE FUNCTIONS A(6d) and B(8) 





tan-!(B(8)/A(5) 


FIGURE 2. 


O Ursell [1949] 
+ This work 


THE FUNCTION 





tan-! (B(8)/A(5) ) 
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N(5,@) = P5686) P57,66,9) + ® na 


t 
iM 


=1 


BI 


cos (2m=+1)0] + ae cd Poe (Ssin®@). 


OUe: 
em=-1 


>. Po bcos 2m + 
m=] 


Here we have substituted the definitions of 6699) and ~ (599) 
from Section 4,2 and a form suitable for numerical calculation of 
~.,6548) from Appendix A-3.4. The open-ended summations over the in- 
dex n in M(5,@) were automatically run to "convergence" on the dig- 
ital computer, where "'convergence'"' means the least significant digit in 


the floating-point calculation. 


The M(65,@) and N(&,0) were calculated for each 6& at inter- 
vals on the cylinder surface never larger than 909 = Oo Cig oo 

The values of M(5,@) and N(6,0) available at @ = 0°(10°)90° 
were used to calculate M (8) and N(8) by a numerical quadrature. 


The defining equations are, from Section 3.8, 


1 t/e 
M65) = af M(5,8) W(@)de 
fe) 


n/e 
f  N(&,8) w(e)ae . 
Oo 


Qik 


NACo) = 


In this case G=l and W(0) = cos@. It is possible in the case of 

the circular cylinder in infinitely deep fluid, as shown by Ursell [1949], 
to substitute the defining equations for M(6,0), N(5,90) and to simplify 
before calculating Mos Nie Here, however, we are concerned with the 
values of M(5,@), N(5,0) which other calculators bypassed. Using the 
calculated values of M(5,@), N(5,@) in a numerical calculation of 


M (6), N¢8) gives an important check since (Section 3.92) 
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The procedure employed was to perform a Lagrangian interpolation 
based on the calculated values of M(5,@), N(6,@) at @= O (ie oo" 
and then to use Simpson’s Rule for approximate integration of 


M(6,@) W(@) and N(6,8@) WC@) with a convenient interval of @. 


The calculations were tested at each 6 by computing 
M(6)AC6) - N (6) B(&) and observing the percentage deviation of 
the calculated value from Taree Curves plotted through error index 
points are shown in Figure 43 as a function of the frequency parameter 
for sample calculations with six and nine expansion coefficients (each) 
Po,6 8) » 45,66) At higher values of 6, seventeen expansion coeffi- 


cients were calculated by the discrete=-angle method. 

Ip, Hydrodynamic Pressure Distribution on the Cylinder Surface. 
The hydrodynamic pressure at the location @ on the cylinder 

surface in units of the hydrostatic fluctuation is (Section 3.7) 


oe , MCS OBS. + N(6 BALE) (unit acceleration) + 
fren NCGS a 


oO 


ae MC5 G2A() =- N(6 9) 5(8) Gind elect) 
A(6)* - B(8)* 


This is a simple, but repetitive, algebraic calculation using the 
principal dependent variables. The pressure fluctuation in phase 


with acceleration, 


MB + NA 
Deaec = > 3 
A” + B 


and the pressure fluctuation in phase with the velocity, 


MA = NB 


Pp, = 6 
Vv Ac Pt 2 
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FIGURE 3. CALCULATION ERROR INDEX vs FREQUENCY. 


2 
(MOA ~ NB) =n /2 
72 
(2) Six expansion coefficients. Method of integration. 


(2) Six expansion coefficients, Discrete angle method. 
(3) Nine expansion coefficients. Discrete angle method, 


Index 2 x 100% 
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were calculated for given 6 and © = O- Clo 3S0.= 


Curves through calculated points of p, (6,8) and p,(65@) are 


shown for example and Tabie I lists calculated values. 


These computed values are used for a prediction of the total 
pressure fluctuation on the cylinder used in the experimental work re-= 


ported in Part C. 


It can be shown by deriving the velocity potential appropriate 
to the case of large 6 that the angular dependence of the hydrody- 


namic pressure appreaches cos®, that is, for large 6 


p(6,@) = 6 cos@ (unit acceleration). 
per, 
Therefore 
a MB + NA 
Go ae) ee = 
A” + B 


The sample calculations illustrate this since it is seen from Figure 6 
PY tends to zero and from Figure 5 that p,/6 tends to cos@ for 
large values of 6. The same result holds for any cylinder of the el- 


liptic family. 


54 The Added=-Mass and Damping Coefficients. 


The ratio of the added mass to the mass of fluid displaced by 


the cylinder (Section 3.91) 


na“p/2 A’ +B 


and the damping-force coefficient (Section 3.92) 


F 7 oT 
maqp/2 Ac + os 
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FIGURE 4. p (8.8) FOR THE CIRCULAR CYLINDER. 
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p (6,0) FOR THE CIRCULAR CYLINDER. 
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p (6,6) FOR THE CIRCULAR CYLINDER. 
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p (6,9) FOR THE CIRCULAR CYLINDER. 
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were previously reported by Ursell [1949] for certain 0 < 6 < 31/2. 
His values are confirmed and the range of frequency parameter is ex- 
tended. The values of M08) N(6) used here were calculated from 
the values of M(8,@), N(6,@) used in the hydrodynamic pressure cal- 
culation as discussed in Section 5.2e2. The results are shown graph- 


ically in Figure 7. 


The results for the damping coefficient can be converted to 


amplitude ratio by 


2 
(eM) 2 OR (2k 
aly. ne 2 a = 


lL 
5 as 658 ?@® «4 


He pointed out in [1957] that although this asymptotic form is correct, 
it may be hazardous to interpolate with this result since earlier theo- 
ries predicted distinct minima. He suggested extended calculations 
could help settle the matter. The present calculations for the ampli- 
tude ratio tend distinctly to the asymptotic value as given by Ursell 


and shown in Figure 8. 





COEFFICIENT 


OS Ursell 1949) 
+ This data 


FIGURE 7. THE ADDED-MASS AND DAMPING 
COEFFICIENTS FOR THE CIRCULAR 
CYLINDER. 
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FIGURE 8. THE AMPLITUDE RATIO A FOR THE CIRCULAR 
CYLINDER. 
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VI. Kliiptic and More General Cylinders in 


Infinitely Deep Fluid 


6.1 Description of the Method. 


Sample calcuiations were made for the hydrodynamic pressure 
distribution, addedemass and damping coefficients of elliptic and 
more general cylinders in infinitely deep fluid. An important ob- 
jective of the sample calculations, in addition to that of obtain- 
ing the hydrodynamic quantities of interest, was to test the flex- 
ibility of the program written for control of the automatic digital 
computer. The program was successfully tested over a wide range of 
cylindrical shapes generally encompassing those of naval architec- 
tural interest. Some of these sample calculations are reported 
here. It is intended that a more comprehensive atlas will be sepa- 


rately reported later. 


The calculations were done on the IBM 704 digital computer at 


the Computer Center, University of California, Berkeley. 
6.11 Numerical Procedures. 


The expansion coefficients Po, 5)» Qo,,66) were calcu- 
lated by the discrete-angle methed described in Section 5.11. The 
expansion equations are given in Section 4.3. The forms of 
Yoq 6519) and Yoga Oot/2) suitable to numerical calculation are 
given in Appendix A=-3.4. The number of expansion equations m and 
the discrete angles e, were varied. The least number m in the 


calculations reperted here was nine so that the hydrodynamic pres~ 


sure could be calculated at ten degree intervals with expansion 
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coefricients known to satisfy the truncated expansion equations 
exactly at ten degree intervals. More often in these calcula- 
tions m = 17. In every case the expansion coefficients were 


tested at arbitrary 6 = 6, # ©, to provide the accuracy index 


t 
@eocribed in Section 5.11. 
The principal dependent variables A(&), B(5), M(6,@), and 
N(6,0) were calculated in the same manner as for the case of the 
circular cylinder but with appropriate definitions of Yama (Oet/2) 
and Po ma 9) as given in Section 4.1. The dependent variables 
M fo), Nf 5) were calculated by numerical quadrature from the cal- 
culated values of M(6,6), N(5,@) as described in Section 5.22 


and the result was tested by comparing MA “ N,8 with ee 


6.12 The Selection of Mapping Parameters 2on.1° 


The number N+ 1 of mapping parameters 
Bone? n=0,1,2,.00e0,N (Sections 2.2, 3.2) has been left unspecified 
in the problem solution. It is convenient here to classify the 
cylinders by the number N+ 1. The naval architect is more ac- 


customed to classifying ship sections by other parameters so some 


of this nomenclature is introduced, 


The ratio of the half-breadth of the cylinder at the free 
surface to the depth of the cylinder at the centerline is designated 
HS 


H = x(n/2) / y(0) = b/d. 


The ratio of the crossesectional area of the cylinder to the area 





oe, 


exinm/2)y(0)v= "2 bd; 


Fas 
Crome oe eran y (en+1)(a, 41) / 204 
Beek au/70 (2nt1)e2 
G n=0 fps 


is designated S and called the "area coefficient." 


The one-parameter cylinders generated by assigning N = 0 
suitable values to ay =O are members of the elliptic family. 


The special case a, = O is the circular cylinder. The ratio 
wee= () + a,) /(1- a,) is uniquely determined. The area coeffi- 
cient for any elliptic cylinder is constant; S = n/4. Thus in 


naval architectural terms the one-parameter family of cylinders 


offers a choice of beamedraft ratio but constant area-coefficients. 


The two-parameter cylinders generated by assigning N= 1 
and suitable values to a, = 0h az =O, to provide more-or-less 
ship-like cross-sections are commonly called Lewis forms, the ratio 
H and the area coefficient S, 


1 + as +a 


Si 1 a “ 
Poe es 
1 - Ree Be 
2 pe a Ee 
Sore Te 2 
+ az - ay 
are determined by ais az and conversely for ship-like forms. Thus 


the two-parameter family of Lewis forms permits some selection of 


area coefficient at fixed beam-draft ratio. 


Still greater flexibility in selecting a geometric form with 


given beam-draft ratio and area coefficient is possible in the 
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three-parameter family. Here N = 2 and appropriate values are 


assigned a L. Landweber and M. Macagno [1959] intro- 


ays az 5° 
duced a systematic treatment of these forms and present a conven- 
ient method for finding the asymptotic value of the added mass 

np (Section 4.4). The introduction of the third parameter specif- 
ically permits a variation in the vertical distribution of the cross- 
section area at given beam-draft ratio and area-coefficient. This, 
in naval-architectural terms, is a two-dimensional interpretation of 


the vertical prismatic coefficient. Ship-like forms with higher area 


coefficient are also possible. 


Prohaska [1947] considered two special cases from the four 
parameter family, N = 3. He selected non-zero values for ass ap 
and for a 


a The examples be illustrates of the case az a 


3° Te 7 


not zero are cusped, while several in the group ass hs, not zero 


may be called ship-like. 


The program for the control of the automatic computations 
permits the selection of any foreseeable number of coefficients. 
Coefficients are stated consecutively up to the highest non-zero 
value and the summations over N terminate:at that number. Many 
sample calculations have been done with various combinations of 


ee@etiicients through 4 The results of sample calculations of 


ar 


examples from the one-, two-, and three-parameter families are pre- 


Sented here, 


6.2 Results for Examples from the One-Parameter (Elliptic) 


Family of Cylinders. 


The sample calculations for the circular cylinder given in 
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Chapter 5 are a special case of the one-parameter family of ellip- 
ero Cylinders The results of calculations for elliptic cylinders 


with H=1/5 and H = 3/2 are presented in the following figures. 


The hydrodynamic pressure fluctuation in phase with the ac- 
celeration, ps and that in phase with the velocity, Py (Sec- 
tion 5.3) are shown for typical values of @ measured in the ref- 


erence plane as a function 6 = Kx(1n/2) = Kbe 


For large 5 >o® =, 


p 
— MC6,0)A(6) = N(5,0)B(4) ne 


hanye a Gs)o 


for all cylinders of this family (Section 5.3). This fact is illus- 
trated by observing the slope of the p,6549) lines and is particu- 


larly evident in Figure 10. 


As pointed out in Section 4,4 the free-surface correction 
kK, = n/m is numerically equal to the ratio of the added mass m _ to 
the mass of the fluid displaced by a circular cylinder of the same 
beam in the case of elliptic cylinders. Therefore curves of k,,(8) 
show the frequency~dependence of the inertia»coefficient of Lewis 


[1929] or the added»mass coefficient of Landweber and Macagno [1957]. 


Curves of the ratio of waveheight far from the @¢ylinder to the 


amplitude of vertical oscillation, 


A = h /h 
are shown in Figure 14. The damping coefficient may be calculated 


from the amplitude ratio (Section 5.4), 
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FIGURE 10. p_ (6,8) FOR THE ELLIPSE H=I/5. 
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Tasai [1959] has previously reported calculations for the free- 
surface correction factor Ky and the amplitude ratio A for the two 
ellipses reported here over a limited range of the frequency param- 
eter. His method of calculating six expansion coefficients satisfy- 
ing eight expansion equations in the least-squares sense was reported 
in Section 5.11. His procedure for evaluating the several integrals 
was not reported. The trends of ali his calculations are confirmed 
and differences within his range of calculated values (determined 
graphically from his figures) are almost certainly due to differences 


in numerical procedures, plus the graphical dispersion. 


These curves give insight into the frequency dependence of the 
various hydrodynamic quantities for forms of various beamedraft ratio 
and constant area-coefficient. The vertical distribution of area as 
measured by the fractional depth of the section centroid or the frac- 
tional depth of the radius of gyration about the free-surface (x-axis) 


m> also constant. 


6.3 Results for Examples From the Two-Parameter (Lewis Form) 


Family of Cylinders. 


The Lewis forms have, for given beam-draft ratio, a range of 
values of area-coefficient. The vertical distribution of area of the 
forms of fixed beamedraft ratio also changes Come, the three-parameter 
family of the next section). This one added mapping parameter fortui-~ 
tously provides sections that are more full than the ellipse and gen- 
erally accepted as ship-like. The upper limit of fullness obtainable 
is somewhat arbitrary. Lewis [1929] preferred to judge the forms by 


the practiced eye of a skilied naval architect; Landweber and 
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Macagno [1957] formalized this decision. The lower limit of fullness 
is probably also arbitrary since the forms tend to be not aestheti- 
cally ship-like and if fine enough, are cusped at B =O for H<¢ 1 


and at BP =n/2 for H Dl. 


The sample calculations made to date suggest that the behavior 
of the various hydrodynamic quantities does not change abruptly as 


the values of ass ee, encompass the reasonably expected range. 


The results of two of the calculations for this family are 
shown graphically. Both forms have H=1, the same as the circular 
cylinder. One form is full and ship-like. The other is the section 


of minimum area and is cusped. 


Figures 15 and 16 show curves through calculated points of 
p, (6,0). The asymptotic behavior of p,/6 for large 6 for this 
family is 
Pa (l+a,) cos 6 - az cos 30 
fe) 


as 6567@. 
(lea, = az) 


This behavior, more complex than simply cos@ as in the one-parameter 


family, is shown by the slopes of the curves in the two figures. 


Figure 19 shows curves of the free-surface correction factor, 
kK), = n/n: It is recalled that for the more general cylinders the 
free-surface correction factor kK) is no longer numerically equal 
to the added-mass or inertia-coefficient. The calculation of Tasai 
[1959] for k, of the full-form is shown over his range 


Oo4+< & < 2.6. 


The damping coefficients for the full form and the cusped form 


are shown in Figure 20. 
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6.4 Results for Examples From a Three-Parameter (Landweber-Macagno) 


Family of Cylinders, 


The use of three parameters a in the mapping func- 


ass az 5 
tion permits some independent variation of the vertical distribution 
of area for the ship-like section of given beam-draft ratio and area 
coefficient. L. Landweber and M. Macagno [1959] chose to measure the 
vertical distribution of area in terms of the parameter n defined 
as the moment of inertia of the actual section about x-axis made non- 
dimensional by the product x(1/2) y (0)? = bd”. Another useful char- 
acteristic of the three-parameter family is that it provides ship- 


like sections of larger area-coefficient than is possible in the two- 


parameter family. 


Examples are shown which illustrate the behavior of the hydro- 
dynamic quantities for a bulb-like section and a section of maximum 
fullness. The sections are shown in Figure 21. The full section is 
an example of sections possible in the three-parameter family that 
are ship-like and have area coefficients larger than is possible in 
Lewis forms. The area coefficient for the full section is S = 1.0ldée. 
This section was selected from Figure 4 of the Landweber and Macagno 
reference cited above. The bulb-like section has an area-coefficient 
of about 0.6 and it is a member of the special group az = O; ays 2 


not zero from which Prohaska [1947] selected examples. 


The amplitude of the pressure fluctuation in phase with the ac- 
celeration at locations near the vertical centerline of the full sec- 
tion increases much more rapidly with increased 6 = Kb than for the 


more rounded elliptic or bulbous forms. 








FIGURE 21. THREE CYLINDERS, H=1/5. 
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The asymptotic form of p,/6 for large 6 for the bulbous 


form is shown in Figure 24, 


In Figure 27, showing the free surface correction factor Ky 
and the amplitude ratio A, data on the elliptic form with the same 
beam-draft ratio, previously given in Section 6.2, are reproduced 
for reference. The amplitude ratio for the bulbous shape indicates 
a relatively high damping coefficient; the full section (larger 
vertical prismatic coefficient) is relatively less damped in heaving 


Ssciilation. 
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FIGURE 24. Pp. (6,0) FOR THE BULBOUS SECTION, H=I/5. 
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FIGURE 26. Pp (6,6) FOR THE BULBOUS SECTION, H=I/5. 
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FIGURE 27. Ky AND A FOR THREE CYLINDERS, H=1/5. 
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Vil. Experiment and Results 


Jel Purpose and Direction. 


Fo Ursell, who obtained the first rigorous sclution for the 
circular cylinder in fluid of infinite depth [1949] also said clearly 


"experiments are needed" and has reiterated this even in recent times: 


"It would be highly desirable to compare these calculations 


with suitable experiments." [1957]. 


The response has not been overwhelming. The earlier experiments 
reported in Chapter I are not strictly comparable, and in later 

years, three-dimensional experiments (Golovato [1957], Gerritsma [1957]) 
have received more attention. Apparently, the only modern experiment 
simulating two-dimensional behavior is that of Yu [1960], who mea- 

sured waveheight ratio and found satisfactory agreement while the prese 
ent experimental apparatus was under construction. Yu reported instru- 
mentation difficuity in measuring the vertical force acting on the 


cylinder. 


Pressure fluctuation measurement was selected as the primary 
objective of the experiment to be discussed. The total pressure 
fluctuation signalled bv a suitable transducer is a distinctive func- 
tion of the location @ on the cylinder, as well as the frequency 6, 
according to the preceding theory. It is considered unlikely that 
reascnable agreement could be expected between experimental results 
and the distinctive prediction unless the analytic model accurately 
represents the essentials of physical reality. Vertical force acting 


on the cylinder was selected as a secondary objective because, first, 
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this 1s an integrated result of a predictable primary physical quan- 
tity. Second, it appears not feasible te avoid encumbering the mea- 
surement of force by ordinary body inertia reaction, which is far 


from negligible, yet not of eoncern to this theory. 


An additionai reason for placing emphasis on the’ pressure 
distribution measurement is the expestation that such investigation 
will help to provide insight into similar problems of naval archi- 


tectural interest. 


The following measurements were obtained with a cylinder of 
10-inch radius in forced vertical harmonic motion of small ampli- 
tude. Details of the experimental arrangement and simulation of two- 
dimensional conditions are recorded in later sections. Over a cer-~ 
tain range of oscillation frequency, 0.2 < f < 0.6 cycles per second, 
the particular arrangement used simulated two-dimensional conditions 
less successfully and this caused consistent small deviations from 
theory in this range. The experiments were conducted at the Ship 


Model Towing Tank, Richmond Field Station, University of California. 


Ls Results of Measurements of the Pressure Fluctuation, 


The total pressure fluctuation at a point on the surface of 
a horizontal circular cylinder given forced vertical oscillations in 
an ideal fluid of infinite depth is predicted with the aid of the 
calculatione reported in Section 5.3. The total pressure fluctuation, 
in units of the hydrostatic fluctuation for a forced motion y=ho sin at, 
is 


C6 0O9t) @ [lep (6,6)] sinws+ p_(6,0) coswt . 
pen. e ' 
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Here p (5,9) and p.\6,0) are the acceleration- and velocity- 
phase components of the hydrodynamic pressure given for this case 


in Section 5.3. The amplitude of the total pressure fluctuation is 


. Ve 
oy BS Ae ee ee 
p(6,@) pen. itilep )” 4 P,, ] F 


The predicted results are shown by solid lines through cal- 
oe 2 7 28 £ ie) e) a ° 
culated points in Figure 25 for @= 0°, 40, and 60. Experi- 
mental points for measurements at the same locations are shown on 
the same figure. These points are (amplitede of the recorded 
pressure fluctuation)/Camplitude of the recorded pressure fluctua- 
tion for hydrostatic displacement}. The distinctive trends, for 
each @, of the predicted and measured results are generally cone 


fiuent. An observed departure from two-dimensional conditions in 


aw 


do 


the experiment is the principal cause for consistent small dev 
tions from predicted results in the range 0.2 < f < 0.5 ecycles/second. 
The trend toward higher-thanepredicted resuites at 6 = 60° at the 
higher frequencies observed may be causec by viscous effects in the 


locations of higher tangential velocity. 


Vo > Results of Measurements of the Vertical Force 


Acting on the Gylincer. 


The tatal vertical force per unit Length required to sustain 


Steady-state motion yz=h. sinwt of the test cylinder is the sum of 


Q 





MB+ NA = 
eee oe re ' 
apa? 25m » apa EE aw 
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and the hycrostatis force caogy and the inertia reaction of the 
mass of the cylinder. The mass of the cylinder in this experiment 
was made equal to that of the displaced fluid for neutral buoyancy 
at rest. The total force in units of the change in hydrostatic 


force is 





Fit) = ‘3 On 
meee = LL} (2+k) 62/4] coswt - == (67/4) sin wt 
2 ¢ 
pear, Mage soe 
, MB+NA 
eens ° 
: ne + ee , 


The amplitude of the total force fluctuation is the square-root of 





(1 -« (Lek) 6n,/4)° + f Ain Convey as ; 
A” + B 


The added mass coefficient k and the damping forse coefficient 


Baa + 3B“) for this case are reported in Section 5.4. 


The predicted result is shown by the solid line through cale- 
culated points in Figure 29. Points for measured values are shown. 
These points are (amplitude of the recorded total force signal.) 
divided by (amplitude of the recorded signal for hydrostatic dis- 
placement). The general agreement is satisfactory. Again, the 
observed departure from twoedimensional conditions is the principal 
cause for smail differences from predicted results in the range 


O.2< £< 0.6 cycles/second. 


Measured results are generally higher than predicted in the 
neighborhood of the force minimum around £ = 0.9 cpa. In this 


region 
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FIGURE 29. CALCULATED AND MEASURED TOTAL 
VERTICAL FORCE. 
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L = (1+k)65n/4 ZO 


and the total force is nearly that in phase with the velocity. Any 
additional energy dissipation (due to viscous effects, end-effects, 
etc.) would contribute additional force increments in phase with the 
velocity and not predicted by the theory. The measured values there- 
fore may well be expected to be higher than predicted and more dis- 
tinctly so when the much larger component in time-quadrature vanishes. 
o4 The Cylinder and The Experiment Arrangement. 

h 


if The Cylinder and The Motion Generator. 
The cylinder was fabricated from a one-quarter inch thick 


plate of aluminum rolled to a circular submerged cross-section with 
10-inch outside radius and le-inch vertical sides above water. The 
cylinder is 41 inches long. The ends are closed by end-plates of 
three-quarter inch aluminum plate. Two heavy steel athwartship 
frames, machined to semicircular shape, together with the end frames, 
assure semi-circular shape,.contribute rigidity, and anchor parts of 
two force dynamometers. The model plus the parts of the dynamometer 
in the model was weighted to neutral buoyancy at the desired water- 
line. The high vertical walls anchor dynamometer parts, contribute 
torsional rigidity and may be used in other experiments with deeper 
Submergence. These high weights cause static instability but the 
overturning moment is small and of little consequence since the model 
is secured by means of vertical-motion dynamometers to the rigid 


motion generator, 


The vertical-motion force dynamometers consist of an inter- 
changeable set of beam "Springs" each 7-1/2 inches long, e-inches 


Wide and 5/32 or 3/l6-inches thick for choice of spring-constant. 
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FIGURE 30. THE MODEL 








FIGURE 31. THE VERTICAL MOTION DYNAMOMETER 
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AS array of aix such aprings connects each of two motiom=generator ree 
ciprocating rods to its steel-frnme amd endoplate im the cylinder in 
am array strongly preferential to vertical relative motion, The rela- 
tive displacement between the model and the reciprocating rod, a meae 
sure of the vertical force acting, was sensed by a Statham Instruments 
Corporation Model G-10-3 umbanded straiuegauge displacememt transducer. 
Relative vertical metions of about 0.002 inch were permitted and this 
is considered a megligible uncertainty in the stroke of the motion 
generator. The auplitude of the stroke was varied from ome-half to 
one imch amc had mo odservable affect om the recorded. measurements, 
Almost all of the preseure measurements were made with eme-half inch 


stroke, 


The motion geaerator is anu elaborate machine desigmed and 
build for the Ship Model Towlag Tank, A 3/4-horsepower consatant- 
speed motor, througe a continucusly-varieble speed redugcer, drives 
& pair of mechanisms patterned after the "scotch-yoke” which transe 
late comustant-apeec rotation to uniform simple-harmomie reciproca- 
tion of twe steel rods. The scotch-yowe, at the expense of elabora- 
tion, has tke advamtage of providing alimost exact sinusoidal motion. 
The motion was continuously monitored by mechasically-drivem poteme 
tiometers on the mechanism providing electrical aiguals which, bee 


cause of orientatiosx, were sine and cosine signals. 


The pressire gauges were inatailed ia the shell of the cylin- 
der at various angular locations at mid-lemgth. The gauges were une 
bonded strain-gauge types with sensitive diaphrame about 3/4-inch 


diameter, The diaphragms therefore subtend (0.75)/10 20.075 radian, 
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about 4.3 degrees and it is understood results measured "at 40 degrees" 
admit the diaphragm extended abcut + 2 degrees. The Statham Instru- 
ments Corporation Type PM-234-TC pressure transducers have calibra- 
tion factors close to 23 millivolts/psi with strain-gauge bridge exci- 
tation of 5 volts. The full-scale range is + 0.5 psi differential. 
The differential feature was used to provide reareside pressurization 
(by air and static manometer) to balance the hydrostatic pressure pre- 
lead on the sensitive diaphragm. This contributes to maximum initial 
sensitivity and linearity; the response of diaphragm-type transducers 
was studied extensively and is reported in the next chapter. The dy- 
namic response of these gauges is uniform from hydrostatic up through 


frequencies well above those reported here. 


742 The Simulation of Two-Dimensional Conditions. 


The experiments were conducted in the Ship Model Towing Tank, 
Richmond Field Station, University of California. The tank is 200 
feet long, 8 feet wide and 6 feet deep. The cylinder and motion 
generator were installed mid-length and transversely between false 
walls installed in the tank to narrow the width at the model to 42 
inches. The model is 41 inches long so there was 1/2 inch clearance 
at the ends. At the frequencies observed, surface action and dis- 
turbances at the ends was not violent but was increasing with fre- 
quency. At frequencies above those reported here the surface wave 


action was not two-dimensional. 


The false walls at the model were 8 feet long, extended above 
the surface, and were solid to tke bottom. The walls were 5-pound 


steel plate rigidly stiffened on the outside by welded flatbar 
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stiffeners. The varailé¢ fa s¢ walls at 


tf} 


4d 


the model were joined to the 


regular tank walis by further sertions at each end about 8 feet long. 


Thus the mcdel was in the 42 insh spare between faise walls which were 


parallel for 8 fee+ and fiared cut on each side to 96 inch spacing. 


It is agreed that for the experiment repcrted here alone, a model 8 


feet long withcut special walis would be preferred, but other work 


was also in mind. Nevertheless, with the exception discussed below, 


it appeared that two-dimensional conditions generally prevailed and 
regular progressive surface waves were generated. Since the motion- 


generator was overepowered and positiveiy driven, steady-state condi- 


tions were approached in the neighborhood of the cylinder very rapidly. 


Under certain conditions, described below, a distinct departure 


from twcedimensicnal canditions was cboseryed. The influence of this 


effect as reflected for example in the total force measurement was 


strongly dependent on the transition conditicns at the ends of the 


parallel walls. Three transition ccnditions are shown in Figure 32. 


at 
iG 


Arrangement (a) is the leass abrupt transition from 42 inch 


Spacing at the model to the 96 insh tank width. This is the arrange- 


ment used for the data recorded in the preseding sections. The cone 


sistent lowerethan-etheoretisal experimental vwaluee in the range about 


0.2 < f£ < 0.6 cycles/eeroni are believed to be caused principally by 


this transiticn, The effest is much mere intense with more abrupt 


transitions. Two obher arrangements were tried and are illustrated 


in the figure tcgetier with ske*ches 
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eliminate this cause of deviation of measured from predicted values. 
In all arrangements the surface wave-action, even in this critical 

range, appeared to be normal and two-dimensional in the near neigh- 
borhood of the model except for minor disturbances in the corners at 


the ends. 


743 Electrical Instrumentation. 


The unbonded strain-gauge displacement transducers in the 
force dynamometer and pressure transducers in the shell of the model 
(Section 7.41) were connected to Brush Instruments Division, Clevite 
Corporation, Model RD-5612-00 carrier and recording amplifiers. The 
carrier frequency is nominally 2 kcps, the bridge excitation Wolters 
normally available is less than maximum usable by the transducers 


but signal strength was sufficient and reliable. 


The electrical signals monitoring the motion (Section 7.41) 
were connected to the directecurrent recording amplifier sections 


of other Brush RD-5612-00 units. 


All electrical signals were recorded on an 8-channel Brush 
Instruments Model RD-2684-50 linear chart recorder. Data were re- 
duced from the records manually; recorded amplitudes measured in the 


range 0.6 to 4 centimeters. 


The frequency of oscillation was computed from the recorded 
motion sinusoids and known chart-speeds and also compared with rev- 
olutions counted by a Berkeley Instruments Corporation "events-per- 
unitetime" counter. This instrument was given electrical pulses 
from cam-operated microswitches and the number of counts in a known 


time interval confirmed the oscillation frequency. 
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VIII. General Study of Instrumentation for Measurement of 
Pressure Fluctuations. 


8.1 General Characteristics of the Pressure Measurement. 


This chapter presents an analysis of pressure transducers 
and their use in measuring the hydrostatic and hydrodynamic pres- 
sure distribution on a body oscillating ina free surface. The 
material is part of a general study of instrumentation pertinent 
to the experimental investigation reported in the preceding chap- 


ter and is of interest in similar investigations. 


The predicted values for the experimental case may be used 
to fix ideas on the pressure fluctuation magnitudes and frequency 
range. The fluctuation is periodic at the model-motion frequency, 


Its amplitude is 


Ww 2 
B(6,0) = ete = ((l-p)* +p], 
(e) 





where P6519), p (5,0) are the acceleration- and velocity-phase 
components calculated for the example in Section 5.3. Curves 


showing this result are given for various © in Figure 28. 


At low frequencies, 65 < 0.2, p(6,0) fale independent of 
location 96. That is, the pressure fluctuation is approximately 
the hydrostatic change due to slow heaving. The surface is barely 
distorted, and the pressure change is nearly independent of loca- 
tion. This is rigorously correct as 6 >0O. Since the displacement 
is assumed to be small with respect to body dimensions, this refer- 


ence pressure change is small. 
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The amplitude of the pressure fluctuation goes through a 

minimum at a characteristic frequency depending on the location 0. 

For example, the minimum is about 20% of the hydrostatic change at 
locations in the bottom sixty-degree sector and occurs at 1.5 < 8 < 3. 
Twenty percent of the expected hydrostatic pressure fluctuation is 

a useful index for the minimum expected pressure fluctuation. For 

the eO-inch-diameter cylinder used in the experimental part of this 
work, ¢O percent of the one-half inch motion in water corresponds 


to about 0.004 psi. 


At higher frequencies, the magnitude of the pressure fluc- 
tuation increases. Ultimately, cavitation at the bottom center- 
line or viscous effects or short surface wavelength violate the 


assumptions preceding these predictions. 


The frequency range of interest in the experiments with a 
cylinder of 20 inch diameter was O-2 cycles per second, that is, 


mgout O< 6 < 3, 


The low frequencies and preference to record hydrostatic 
pressure suggest a diaphragm-type transducer. Piezoelectric-type 
hydrophones are contra-indicated by this requirement because of 


their limited low frequency response. 


The low pressures and pressure fluctuations dictate high 
Sensitivity. For ship-like cylinders of about 1-foot beam the 
minimum detectable pressure change should be about 0.001 psig (about 
1/40 inch water or +110 DB ref 0.0002 AS/en ee This requirement 
is easily met by a piezoelectric transducer but the high sensitivity 
conflicts with the requirement for small diameter in the case of 


diaphragm types. 
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For experimental ccnvenience, cylinder dimensions must be 
manageable and for maximum utility the transducers should be com- 


patible with the ship models ordinarily used. 


8.2 Piezoelectric Transducers. 
8.21 Open Circuit Operations the Piezcelectric Wafer. 


The piezoelectric transducer gives promise for this applica- 
tion because of its sensitivity. Its disadvantages are limited low- 


frequency response and sensitivity to local vibrations of its mount. 


To illustrate the cperation of a piezoelectric transducer 
and its major disadvantage, consider the response of a disc to nor- 
mally incident pressure signals. The disc area is A, its thickness 
t. Assuming a material similar to nalaen titanate, the disc is po- 
larized across the thickness direction, a. The rear face is blocked 


by a relatively large mass. 


This wafer of piezoelectric crystal, compressed in the direc- 
tion of its thickness + with a pressure p toa stress level 
p= F/A, generates a proportional electric field in the direction 
of polarization: | 


E KpR , 


ii 
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The potential difference across the water is 


Vrs (kt)o ="op 


7 aie 


G 
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pressure sensitivity, volts/microbar. 





oh) 


The constant k in units of volts/cm thickness/microbar pressure 


is typically 1.41 x 107? 


2, 


for ceramic "B'' Barium Titanate and 
2.41 x 10 ~ for common piezoelectric materials called PzT-4 ana 


Bate5. (1 microbar = 1 dyne/em“). 


The open circuit sensitivity of a piezoelectric transducer 
is therefore independent of diameter and is easily of the order of 
-110 db referred to 1 volt per microbar., Thus piezoelectric trans- 
ducers easily meet the required characteristics of high-sensitivity 


and small diameter. 


The required low frequency response requires closer examina- 
tion due to the loading effect on the crystal of the external cir- 
cuiltry. To examine this, cne uses the equivalent circuit of the 
crystal element. This consists of a voltage generator of zero in- 
ternal impedance in series with an impedance equal to the internal 
impedance of the crystal. The voltage output of the conceptual 
generator is equivalent to the crystal pressure sensitivity. The 
internal impedance of the crystal is considered constant and equal 
to that of the capacitor formed by the crystal of dielectric mate- 


rial between equal-potential conducting surfaces. 


erie 
ec) 
e(t) = op(t) - 


Smee Low Frequency Response of a Piezoelectric Transducer 


and Amplifier. 


The input impedance of an amplifier connected to a piezo- 


electric transducer degrades the low-frequency response relative 
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to the open-circuit performance. This is seen by considering the 


amplifier input signal e, 3 


e(t) 


In the usual circuit-analysis technique we use 


s = jw 
E.(s) = transform of e,(t) 


E(s) = transform of e(t) 


= o& P(s) 
_ Soe £2 SE 2h 
ECs) = E(s) s+ 1/t o P(s) s+ l/r 
oo om 81S oa 


Consider a step-echange in pressure: 


SerOrt < © 
pGt) ={s t 20 és 
The amplifier input signal is 
e.(t) = o ae */t (6.2 0)e 


The open-circuit (unloaded) response of the crystal transducer 
would have been o A. The penalty imposed by the external ampli- 
fier (load) is the exponential decay of the hydrostatic response 


with a time constant +t. 
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e.(t). Time constant 
al 
T= eRe 


Response to a hydrostatic change in pressure, 


The steady-state response to a sinusoidal time-varying pres- 
sure decreases in amplitude and is increasing in phase error below 


mae cutoff frequency. The cutoff frequency is 


Ww = 1/r = 1/Re ° 


Cc 
O 
-10 
Phase error, 
-20 degrees lag 
Relative 
Response, 
db -30 
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Amplitude Response and Phase bkrror relative 
to Open-Circuit Response of a Piezoelectric 
Crystal for Sinusoidal Pressure Fluctuation 
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FRACTIONAL AmptituDE RESPONSE 
OF A PIEZOELECTRIC TRANSDUCER AND 
DC. AMPLIFIER AS A FUNCTION OF 
FREQUENCY NEAR CUTOFF Wet 


FIGURE 33. 
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These resuits show that the itcw frequency response of the 
piezoelectric transducer plus amplifier is limited by the amplifier 
input impedance R and the crystal equivalent capacitance c. The 
amplifier input impedance is bounded by practicai considerations of 
the associated thermal noise and by susceptibility to extraneous re- 
sponse due to the high impedance level. The crystal equivalent ca- 
pacitance c can be increased by increasing the wafer diameter, 
which conflicts with small size; or by decreasing wafer thickness, 
which decreases sensitivity; or by stacking thin wafers, which in- 


creases complexity of the transducer. 


Therefore it can be expected the use of piezoelectric trans- 
ducers will be met with increasing difficulty as the desired low 


frequency response is extended. 


Some techniques are suggested to partly alleviate this limi- 
tation. 
8.23 Extending Low Frequency Response by Capacitive Loading. 


The low frequency response of piezoelectric transducer plus 
amplfier may be extended at the expense of sensitivity. The ampli- 


fier input resistance R is shunted with a capacitance C, 
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Comparing this result with that of the amplifier alone, 


(i) The apparent sensitivity is reduced by a factor 


c/(C+c) < 1.0 


(1) The hydrostatic time-decay factor is increased by 
a factor (C+c)/c. This is equivalent to reducing 
the low-frequency cuteff frequency by the factor 


c/(C+e) < 1.50 


With this technique the relative response is kept uniform 
to a new lower cutoff-frequency at the expense of being at a lower 
absolute level. The degree to which this can be pursued thus de- 
pends on the original crystal sensitivity and the magnitude of 
pressure fluctuations which are to produce measureable signals withe- 
in the absolute level capability of the available amplifier. That 
is, the signal to noise ratio at the reduced apparent sensitivity 
will limit an attempt to indefinitely extend the lowefrequency re- 


sponse by this technique. 


8.24% Extending Low-Frequency Response by Feedback 
Amplifier Technique. 

Manipulations of the low-frequency cutoff and of the appar- 
ent sensitivity, interdependent in the preceding technique, can be 
largely separated at the expense of increased electronic complexity. 
An operational amplifier of standard design is introduced with ca- 
pacitive feedback, The amplifier will have a large gain, -A, with 
negative sense meaning phase reversal. The amplifier is realisti- 


cally assumed to have constant gain over frequencies of interest 
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and to have negligible output impedance. The input impedance cannot 
be neglected and as before is R. This amplifier is to be considered 
aS a preamplifier or buffer for operational purposes between the pie- 


zoelectric transducer and the usual recording system. 





By the usual circuit-analysis methods and without further 


assumptions, 
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Comparing this result with those preceding it is seen the sensitivity 
and time constant are changed by functions of C and A. If the 
feedback element C is reduced to zero the relations are simply 
those for an amplifier following an amplifier-loaded piezoelectric 


transducer. If C is not zero and the realistic assumption is 





made that 
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[A] >>1 


then the relations are simplified. Compared with the amplifier- 


loaded case, 


a) 


(Gr) 


the sensitivity is modified by the factor c/C 


the time constant is increased by the factor 


COLA) c= CGAY oc: 


This is equivalent to reducing the low-frequency 


cutoff by a factor c/CA. 


{It is striking to observe then that if C is made approximately 


equal to 


(i) 


Ga) 


Cy 


then 
the sensitivity is equal to the open circuit 
Sensitivity .« 


the lowefrequency cutoff is reduced by the factor 


-1/A, where A is a large number. 


It appears this should be a useful technique for extending low- 


frequency response with nearly independent control of the appar- 


ear semsitivity and the extent of cutoff reduction. Practical 


limitation on the extent of this manipulation is not clear. 


Be 5 The Pressure Sensitivity of Diaphragm-Type Transducers. 


8.31 


Definition of Sensitivity. 


The transducers considered consist of two elements. 


ED, 


the mechanical pressure receiver. Primarily this 


will be a diaphragm but may be a tube, bellows or 
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FRACTIONAL RESPONSE. AS A FUNCTION 
OF (ADDED / INTERNAL) CAPACITANCE 


FIGURE 34. 
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capsule which deflects under pressure and seals the 


device against fluid ingress. 


ad) the deflection detection schema. This converts the 
pressure-recsiver (diaphragm) deflection to a signal 
in useful form. In order to take advantage of elec- 
tronic (linear) amplification, electrical schemes 
are of most interest in devices of high sensitivity 
and resolution. Examples may be found using various 
physical principles: variable resistance, capacitance, 


reluctance, light Tluxe ete, 


Pressure variations on the pressure receivers of several different 
types of transducers change the internal energy in the detection 
device a negligible amount compared with the change in elastic strain 
energy in the pressure receiver. Limiting consideration to such 
stiffness-controlled transducers, it is possible to analyze the two 
elements of the transducer separately, and the transducer pressure 


sensitivity is the product of the sensitivities of each element: 


5 = 5(p), the diaphragm deflection at the center; 
2.5 e (5), the transducer signal output; 
de. dé de. 





The effectiveness of the pressure measurement system cannot be 

judged on transducer sensitivity de ./ap alone because this signal 
Will be contaminated by a noise signai. Also, a response index will 
be required characterizing the least readable amplitude fluctuation 


of the indicating needie or shart line. Before attacking these 
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problems it is fundamental to study the deflection-pressure rela- 
tion 6 = 5(p). This section will consider this relation and the 


pressure sensitivity, 


16 


Go = ap 0 


8,32 Small Deflections of Clamped-Edge Flat Circular 
Plates; Kircnnofr Theory. 
The theoretical maximum pressure sensitivity of a diaphragm 
is generally considered to be the initial sensitivity of a flat 
Circular plate under uniform pressure load on one side producing 


infinitesimal inextensional strain. 


By inextensional strain is meant that the plate neutral 
plane is not extended. The internal strain energy is greater and 
the plate is stiffer if extensional strains exist. Inextensional 


conditions are consistent with Kirchhoff plate theory. 


The common and useful practice of diaphragm corrugations 
invariably introduces iocations of extensional strain. Corrugated 
diaphragms realize less semsitivity than the theoretical initial 


sensitivity of the equal flat dise in Kirchhoff plate theory. 


The Kirchhoff theory, neglecting stretch and shear of the 
plate neutral plane, yields the following result for the deflec- 


tion 6& at the center of a clamped-edge flat disc of diameter D 


uniformly loaded by a pressure p of one sidé: 


where 
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2 
F = a plate modulus E/(1 -p AR aWoe Bt /12(1 - p) ' 


rs 
iW 


thickness, 


The deflection is linear in pressure; the pressure sensi- 


tivity is independent of the deflection, 


Non-dimensional expressions for deflection and for pres- 


sure are: 
kk 
G . Pa By (2 


These relations are illustrated graphically in Figure 35 
and 36. These representations are subject to the assumptions of 
their theory. Thus it is clearly not correct that 6/t increases 
linearly indefinitely and it is not possible to indefinitely in- 


crease sensitivity oF/D by decreasing t/D. 


We will consider several cases illustrating that the pres- 
sure sensitivity actually achieved by a thin diaphragm will prob- 


ably be less than this theoretical maximun. 


8.33 Pressure Sensitivity of Flat Circular Plates with 
Deflection not Infinitesimal. 

For practical purposes it is not always possible to limit 
diaphragm deflections to deflections small in comparison with dia- 
phragm thickness. Such factors as hydrostatic preload or practical 
deflection detection schemes lead to larger deflections and, as we 


will see, possibly decreased sensitivity. 


The exact solution of the equations of equilibrium of a flat 
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NONDIMENSIONAL DEFLECTION OF 

OF FLAT CIRCULAR DIAPHRAGMS 

WITH CLAMPED EDGES AND NO 
INITIAL. TENSION 
(1) KIRCHHOFF 


(2) WAY 
(3) HENCKY 


FIGURE 35. 
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NON-DIMENSIONAL PRESSURE SENSITIVITY 
OF A FLAT CIRCULAR PLATE 


FIGURE 46 
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circular disc uniformly loaded on one side and with clamped edges 
and no initial tension is due to S. Way [1934]. This solution is 
available in tabular or graphic form but not in algebraic or series 


form. This solution is also illustrated in Figure 35. 


Approximate solutions are due to A. Nadai [1925], 
S. Timoshenko [1940], and A. Griffith [1927]. Another approximate 
solution is reported by E. Waters in the discussion accompanying 
the solution by Way [1934]. These approximate solutions are ex- 
pressed in algebraic form. The deflection-pressure function is non- 


linear and additional definitions of sensitivity are useful. 


Pressure sensitivity is d6/dp as before and may be called 
the instantaneous pressure sensitivity. Initial sensitivity is de- 
fined as the limit of the pressure sensitivity for infinitesimal 


deflections: 





The ideal flatedisc case has constant pressure sensitivity 
d6/dp equal to the initial sensitivity 4d6/dp] and average sen- 
p=0 
sitivity b(p)/P, « 
The graphical presentation of the large deflection case sug- 


gests that the initial sensitivity is equal to that in the ideal case. 


The pressure sensitivity decreases for increasing deflection. The 
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average sensitivity at a given pressure, from algebraic expressions 
Similar to those of the approximate solutions of Timoshenko, Nadai, 


Griffith or Waters is 


oF (pF rc 1 
(2)? 0.545 (amma) (=) 
D D 
[the constant 0.545 may vary about + 10 percent depending on Poisson’s 
ratio and the approximation, Timoshenko reports 0.583 for P= 0.25 
and 0.488 in another approximation. Our constant happens to be an 
average and agrees with Griffith for p = 0.29]. This is illustrated 


in Figure 36 for two assumed valves of 
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It is seen the average sensitivity may weil be markedly lower than the 
theoretical maximum. The instantaneous pressure sensitivity is always 


less than the average sensitivity at finite deflections. 


Corrugated diaphragms such as those of the National Bureau of 
Standards series realize average sensitivities somewhere between the 


ideal maximum and the value of the average sensitivity for 6(p)D = 17 100n 


This limitation of sensitivity to values less than maximum 
theoretical applies when deflections are approaching O.4t or more. 
High sensitivity requires small <+/D ratio for a given material. 
Measurable deflections which are a small part of t force larger D. 
Fortunately, sensitivity is proportional to the fourth power of diam- 


eter. Conversely, attempts to reduce diameter are penalized by this 
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power. An attempt to achieve high sensitivity at given limited 
diaphragm diameter by reducing diaphragm thickness t (i.e., small 
t/D) will ultimately face this limitation when the detectable de- 
flection is of the order of t. Another problem, described in the 
next section, is a loss in sensitivity characteristic of thin dia- 


phragms. 


8.34 The Large Deflection of a Thin Diaphragm with 
Zero Initial Tension; the Solution of Hencky. 

The following comparisons of thin diaphragms, or membranes, 
with circular plates having bending stiffness will use the center 
deflection as a basic measure. It is well to observe that this 
may be an incomplete criterion, depending on the means employed to 
detect diaphragm distortion. Center deflection is sufficient if 
the detection system depends only on this, or if the deflection 
form of the cases compared is similar. The deflection form of stiff 
plates and membranes are not geometrically similar. Therefore the 
volume sensitivity, the volume swept by the deflecting diaphragm 
per unit pressure, is different in each case. A practical case in 
which this would be significant is when the change in capacitance 
is used to detect deflection. A membrane, sweeping out a larger 
volume than a stiff plate at the same center deflection, produces 
significantly greater unit change in parallel-plate capacity (see 
for example Lilly, Legallis, and Cherry [1947]). Since the primary 


index is center deflection, attention will be given only this item. 


One may suspect that a sufficiently thin diaphragm would have 


no effective bending stiffness, so that deflections would develop 
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only membrane tension stresses. Kirchhoff theory and the solution 
of Way do not apply. H. Hencky [1915] has achieved an exact solu- 
tion by power series methods for the equations of equilibrium of a 
flat disc uniformly loaded on one side, subject to the assumptions 


ae mana No initial tension. 


A commonly used plate modulus defined as 


cae oe 
12 PAG ET a 


is in effect reduced to negligible value by Hencky’s assumption. 


Hencky’s result for deflection at the center of this disc is 


won Gbe : eee 


cry m 


(The constant 0.662 is given by Hencky for Poisson ratio p= 0633 


other authors have values + 5 percent or less]. 
This result is essentially predicted by the approximate solu- 
tions for the large-deflection case discussed above with the added 


assumption 


5(p 
t.3 O t 
ee ONS 


so that from those approximate solutions 
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or in equivalent form 
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> = 0.68 poe, 
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Hencky’s solution has the elegance of exactness; this coincidence 


is given as a suggestion of reasonable continuity between solutions 


recorded here, 


Hencky’s solution is also shown graphically with the pre- 


ceding solutions in Figure 35. 


The instantaneous sensitivity is conveniently found from the 


Hencky solution by logarithmic differentiation, 


dé _ 1 dp 
eo 3°” 
Bo chat D 
0.662 i 
oe ft -LS_LL Et _16 
Go, 3p 3 p ° 


The initial sensitivity is indicated to be indefinitely large, 
which is apparent from the slope of the deflection curve. Unfortu- 
nately the very high initial sensitivities appear to be not realizable. 
A physical interpretation of the infinite initial slope of the deflec- 
tion curve may be that the no-load rest position of a thin diaphragm 
without initial tension is unstable. This is possibly a part of the 
phenomena called "oil-canning.'"' Further, if one draws on the analo- 


gous result from the approximate solution, the necessary condition 


2 
&(p 


3 ) 
ie O t 





is not fulfilled for practical (t/D) at sufficiently low 6(p.). 


The assumption for this case implies large (D/t), small 
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t/D, so that very small pressures are equivalent to values of 

pp /et* producing operating points of 6, p away from the ini- 

tial region. A small hydrostatic preload does this for instance. 

For finite deflections the instantaneous sensitivity decreases rapid- 
ly to values less than the constant value of the ideal disc or the 
corresponding values of the large deflection (Way) solution. The 
theoretical average sensitivity of the Hencky-type diaphragms would 
appear to exceed that of the Way-type diaphragm. This is considered 
illusory because the instability of the no-load initial position of 
Hencky diaphragms would in practice mean an average sensitivity mea- 


sured from some initial non-zero deflection. 


This section has employed the Hencky solution to show that 
exploiting thin diaphragms without initial tension and in the limit 
having entirely negligible bending stiffness is not likely to lead 
to practical sensitivities equal to the theoretical sensitivity of 
the Kirchhoff diaphragm. A characteristic of this limitation of 
sensitivity is the assumption of no initial tension. We consider 


now whether there is benefit in adding initial tension. 


8.35 The Small Deflection of a Thin Membrane with Initial 
Tension. 

An exact solution to the equilibrium equations for small de- 
flections of a thin circular disc uniformly loaded on one side and 
subject to a large initial tension T is known. The formulation of 
this problem assumes that the bending stiffness of the plate 
(F = Beo/ie( en) is negligible. It is also assumed that the 


initial radial stress (initial tension) is sufficiently large so that 
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the pressure load is brought to equilibrium by opposing tension com- 
ponents due to the curvature of the membrane and not by additional 
stress. This solution is that of the membrane equation 


7 ~8(r) = 2 i 


which yields for deflection at the center 
2 
D 
5 = ae ° 
The deflection is linear in pressure; the pressure sensitivity is 


constant. 
In order to compare this with previous results, we make sub- 
stitutions, with ES equal to the initial unit strain, 


2 


Do E T-p 
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This is the same form as that for the ideal disc times the bracketed 
factor. That factor, for practical values of (t/D), reduces to 


unity for initial strains e. of the order ost or less. 


The pressure sensitivity of the membrane is constant 


es 
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If one compares this form with the result from Hencky’s solution, 

the factor 3(1 + pie /16 and 0.545 (a7) = are seen to have equiv- 
alent roles in reducing the non-dimensional pressure sensitivity 
oF/D. The apparent high sensitivities of the Hencky solution at 
small 6 are not realizable nor are those of the membrane solution 
for small initial strain since in either case the assumption inherent 


in the appropriate exact solution is not satisfied. 


We will now consider that the Hencky solution may represent 
one extreme condition of a diaphragm with no bending stiffness and 
the membrane solution another. The Hencky solution is for zero ini- 
tial strain. The membrane solution is for such large initial strain 
that radial stress may be considered constant. freee on between 
should be considered since both Hencky and Membrane Solutions give 
the illusion of initial sensitivities greater than that of an ideal 


Kirchhoff diaphragm. 


8.36 The Large Deflection of a Thin Diaphragm with 


Initial Tension. 


The membrane solution is not valid for small initial tension 
and Hencky’s solution includes no initial tension as a boundary con- 
dition. An exact solution of the equations of equilibrium fora 
flat circular disc, assumed sufficiently thin so that nee ti or 
in effect so that the bending stiffness is negligible, and subject 
to initial tension or initial strain, is reported in Appendix B. 
This solution bridges the gap between that of Hencky and the mem- 


brane solution. 


As a result of this solution 
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Here K for a given plate is a function of initial strain and the 
pressure (that is, the finite deflection) and is shown graphically 


in Figure 57. For zero initial strain 


K(O) = 0.662; 


that is, Hencky’s solution is reconstituted. K(e,) decreases in 
value for increasing Eos For large initial strain one has the 


asymptotic solution 
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The membrane solution is asymptotically recovered as it should 
be for large values of initial tension. 
The average sensitivity is 
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For initial tension zero the term K? /p tends asymptotically to 


3 B, 
K 0.662 

5 5 as p 0 
Delle 2 


but this is not useful on the basis of previous remarks because of 
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small residual bending stiffness or because of initial tension 
essential for stability of very thin plates. In case of initial 


tension T not zero, the constant tends asymptotically to 


2 4 p 9 
p T3éO ee Hane 


recovering the membrane action. 


8.37 Composite Approximate Solution for the Finite 
Deflection of a Clamped Edge Diaphragm under 
Uniform Load, 
It remains to state a composite approximate solution for the 
finite deflection of a diaphragm taking into account bending stiff- 


ness, the influence of finite deflections and admitting initial ten- 


sion. 


Griffith [1927] reports carrying out an iterative procedure 
starting from the Kirchhoff (exact) solution for infinitesimal de- 
flections in the general equilibrium equations. His result in our 


algebraic form is 


OF 3 ai 

D 25 S ce . 
t t. (14+ n) (173 -73p) -5 ta Slee) 
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An argument will be made that this form, adjusted, is a correct one 


for an empirical correlation of the sensitivity of general diaphragms. 


The form reduces to the exact (Kirchhoff) solution in the ini- 


tial case with the deflection 6 and initial strain Eo made zero. 


For initial strain E made zero and neglecting the cube of 





124 


t/D this form is comparable with the exact solution of Hencky and 
can be made coincident by small adjustment of the constant 


+ Pp) (173 - 73 pr.)/360. 


With only the initial strain term neglected this form is con- 
Sistent with approximate solutions of the large deflection case 
which in turn are useful algebraic forms representing the exact so-~ 
lution of Way. The constant (1 + Pp»? (173 - 73 p)/360 may require 
adjustment but not to an extent practically conflicting with the 


preceding adjustment. 


For large values of initial strain and neglecting the dis- 
placement and (t/D)? terms this form does not reduce to the 
exact membrane solution without a ten percent change in the constant 
Bil + psa to 45(1 + p)/240 = Clee pd/16 and becomes exact with 


this change. 


Neglecting only the term (t/D)? representing bending stiff- 
ness, the sensitivity is controlled by finite displacement (Hencky) 


or initial tension (membrane action). 


The particular advantage of this empirical form is that it 
makes clear that the pressure sensitivity of a flat disc diaphragm 
will not be greater than the theoretical maximum of the Kirchhoff 
theory and displays that, for very thin diaphragms, the loss of 
sensitivity due to practically unavoidable values of 6/D or E. 


will be considerable. 


8.4 Dynamic Response of Flat Circular Diaphragms. 


8.41 Error Due to Dynamic Response Characteristics of 


the Diaphragm. 


The preceding investigation of the pressure sensitivity of 
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circular diaphragms was limited to static pressure loads and did not 
consider dynamic effects. The dynamic response of the transducer 
Should be considered if the pressure fluctuations include components 
at frequencies comparable with some fraction of the lowest signifi- 


cant frequency in the response of the diaphragm. 


A convenient model for estimating dynamic error is to consider 
the pressure receiver to be a one-degree-of-freedom linear mechanical 
system with a natural frequency equal to the lowest natural mode fre- 
quency of the diaphragm in free vibration. With this model the mag- 
nitude of dynamic error at pressure fluctuation frequencies oe small 
in comparison with the diaphragm lowest-made frequency WO, is pro- 
portional to (w/a) and is minimum for a damping constant about 


70 percent of critical damping. 


For frequencies encountered in the experimental part of this 
work it is not difficult to obtain ratios W/W, that are small. 
For other model sizes or more violent oscillations this may be more 


difficult. 


The high natural frequency of clamped stiff plates or stretched 


diaphragms is helpful in reducing this error. 


8.42 The Lowest Natural Mode Frequency of Circular Diaphragms. 


The dynamic response in vacuum of a circular flat diaphragm 
with clamped edges is studied in W. Mason [1942], The solution is 
applicable to diaphragms with given values of initial tension. The 
result is not conveniently expressed algebraically. A convenient 


graphical solution for the lowest mode frequency is given by 
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2 f 2 
J. Patterson [1952] in the form Ww er d pi -p eit VEg as func- 
‘ 2 a 2 
tion of lee r (l-po)/t > This function is shown graphically in 


Figure 38. 


The least value of the lowest mode frequency for a given dia- 
phragm occurs for zero initial strain. For this condition the above 


form gives 
w are 15 1 -p°) /tVEg = 10.22 . 


This corresponds to the clamped edge stiff plate equation for the 


lowest mode frequency found in handbooks as 


where V is the speed of sound in the diaphragm material. [The 
constants by various authors may vary less than five percent]. 

This form will probably be more convenient than the preceding be- 
cause of the practical difficulty of knowing E and, particularly, 


Ee 
.) 


A simplified form is also available for the case of large 


initial tension, corresponding to the membrane solution, 


w r/V = 2,4 VT/Et . 


The assumptions essential to the latter, or the complete 
solution first given, show that this value of wo, must exceed that 
of the zero initial tension case. Sample calculations of the 
clamped-edge stiff-plate case show the conveniently high values of 


Ww, for small gauges if one takes the velocity of sound in thin 





2Rw, V3p(I-2) /tVEg 


te7 





i2 R%.( j-p*)/ 1° 


GRAPHICAL SOLUTION FOR THE FUNDAMENTAL NATURAL 
FREQUENCY OF A CLAMPED, CIRCULAR, FLAT DIAPHRAGM 
(FOLLOWING PATTERSON [1952] ). 


FIGURE 38. 
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A a 


2 
plates (Eg/pCl = p ) to be for example about 17,000 feet per 


second in steel, 


8.5 Pressure Gauges Mounted Internally. 
8.51 Advantage and Disadvantage. 


The pressure sensitivity of a diaphragm-type pressure re- 
ceiver is a strong function of the diameter, The required pressure 
sensitivity 1s determined by the pressure fluctuation amplitude that 
is expected and the minimum deflection that is detectable. Further, 
the required diameter may exceed the maximum physical dimension al- 


lowable on the surface of the oscillating model. 


The apparent physical diameter on the body surface can be 
reduced by mounting the pressure gauge diaphragm internally. There 
is necessarily then an internal instrument cavity and connecting 
passage or tube. These introduce two effects which could require 
compensation in the measurement; they are (i) the effect of accel- 
eration of the enclosed fluid on the indicated pressure; (ii) at- 


tenuation of pressure fluctuations in the tube and cavity. 


8.52 The Effect of Acceleration of the Enclosed Fluid 
on the Indicated Pressure. 
The fluid enclosed in the instrument cavity and tube is forced 
to accelerate with the body. Therefore there is a pressure gradient 
in the enclosed fluid due to body acceleration. The top surface of a 


cavity is a suitable location for a horizontal diaphragm. 


As an example, consider an instrument cavity in the vertical 


flat side of a vertically oscillating wall. The fluid is assumed 
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AN INSTRUMENT CAVITY IN AN 
OSCILLATING WALL. 


FIGURE 39. 
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inviscid, incompressible, and at rest. Because the wall and motion 
are vertical the fluid external of the instrument cavity remains un- 
disturbed. The external hydrostatic pressure is instantly communi- 
cated to the top of the instrument cavity by an infinitesimal port. 

The external hydrostatic gauge pressure at this port is Po + peldy(t)IJ, 
where 6y(t) is the oscillatory motion and Pp, is the hydrostatic 


pressure at the initial position. 


There is a pressure gradient in the fluid in the instrument 
cavity due to body acceleration a = Heyes « The pressure in the 


cavity is 
p(n) = ple-a)n + Py, + pE FeyCt a 


where n is the vertical distance down from the cavity top to the 


horizontal plane at which the pressure is Po(n) - 


For a top port the indicated pressure on the top surface is 


a direct measure of the external pressure. 


8,53 The Dynamic Response of the Internally-Mounted Gauge. 


The question arises: what is the dynamic pressure sensitivity 
of the pressure receiver to pressure fluctuations at the entrance of 
the connecting tube? It is clear that the hydrostatic pressure sen- 
sitivity is not affected, but it is also well-known that dynamic ef- 


fects may appear. 


The analytic model considered here will exploit the conclusion 
found in Section 8A that the lowest natural frequency of a flat cir- 
cular diaphragm will probably be high with respect to frequencies of 


interest in experiments visualized in this type work. We will now 
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assume that the lowest natural frequency of the diaphragm is also high 
with respect to the acoustic effects of the tube, cavity and resilient 
diaphragm. This assumption enables us to "uncouple" the dynamic re- 
sponse of the diaphragm, which accounts for diaphragm mass and elas- 
ticity, from the response of tube, cavity, and elastic (mass-less) 
diaphragm. Thus we expect dynamic behavior of the tube, cavity and 
elastic diaphragm to be of interest at frequencies low enough to con- 


sider the diaphragm to be stiffness-controlled. 


The problem is to consider the attenuation of oscillatory pres- 
Sure fluctuations in a viscous fluid traveling through a connecting 
tube to an instrument cavity terminated with an elastic but mass-less 
diaphragm. This analysis is reported in Appendix C. Sinusoidal gauge- 
pressure fluctuations of frequency w are assumed at the entry port 
of the connecting tube of radius r. The pressure fluctuation at the 
diaphragm is assumed to be equal to the pressure at the cavity-end of 
the tube, The diaphragm elasticity and cavity volume provide a condi- 
tion on fluid continuity. The result is presented as a complex ratio 
of the indicated- to external-pressure fluctuation P,/P, » The magni- 
tude of the ratio is the amplitude response normalized with respect to 
hydrostatic response, and the argument of the ratio is the phase lag 
of diaphragm sinusoidal response. 

The response P;/P, is characterized by the value of the dimen- 
Sionless parameter war /d . This parameter also characterizes the ve- 
locity profile in the tube. (See L. Loitsyanskii [1957], Section 86, 


H. Schlichting [1955], Chap. XI.) 


For values of is less than one, viscous forces are predominant 
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FRACTIONAL RESPONSE AND 
PHASE LAG OF INDICATED PRESSURE, 
VISCOUS REGIME.. 


FIGURE 40. 





133 


and inertial forces negligible. In this case 


se at eats 
Py 1+ jwrt 
where ne Bia? 
7 n 
2 aaets 
wo =a characteristic frequency constant fora 

given fluid, tube, and cavity, and diaphragm 
elasticity. 


This result is also derived directly with a model assuming uniform 

Poiseuille flow throughout the tube and neglecting inertia effects. 

The relative response is down 29 percent at thecutoff frequency 

w, =1/r and the phase error is 45 degrees lag. The cutoff fre- 

quency is decreased with increased diaphragm elasticity (sensitivity). 
For values of wae /v of about four inertial effects are sig- 


nificant. A convenient approximation is to present the system re- 


sponse in the form 


Ps 1 


oe RGU ae. oo 


Here u is the normalized frequency ratio w/w,, where w, isa 
characteristic frequency analogous to the resonant frequency of a 
moderately damped spring-mass system. For W= Wo the indicated 


pressure lags the driving source 90 degrees. The damping ratio is ¢. 


For large values of a the response ratio tends to 


Ps iL 


eee = ee See e 


Po 1- (w/w,)* 


This is useful in showing the undamped (non-viscous) nature of the 
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acoustic response and also alerts one to the fact that strong acous- 


tic resonance is predicted which might damage a sensitive diaphragm. 


The most general expression is 


4 o) 


se aL 


oh " ~al, Coen al + bL sinh aL 


where 


aL = j 2n(L/A) V1/01 -J) 

bL = ~ (w/w,)*/(1 - J) 

l= 2J, (ka)/(ka) J, (ka) 
(ka) = jua@/y 


For tube lengths L small with respect to the acoustic wave- 


length in the tube A, this tends to 





and approximations for 1/(l-J) give the reductions cited above. 
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APPENDIX A 


Mathematical Notes Referenced in the Text 


Awl Asymptotic Behavior of 95,(%o¥ eh) for 


| x | > OO, Y = O- 
A-2 Asymptotic Behavior of ~(K,h) for 


|x | >@, y = 0.2 


A=3 Expanded Forms Useful in Numerical Calculation 


of Certain Definite Integrals. 


A-3.1 The Coefficients of P56 5) + Go 8) in 


Expansion Equations for the Method of Integration. 
A-3.2 Two Special Forms. 


A-3.3 Expressions in the Right-Hand Side of the 
Po,,6 6) Expansion Equations. 


A-3.4 Alternate Forms for 9. (5,0), Yoq6 519) 


K=S55 Expressions in the Right-Hand Side of the 


d5,(8) Expansion Equations. 
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APPENDIX A 


A-l Asymptotic Behavior of P56 XV 9h) POY || =a ey Oe 


The asymptotic behavior of the function 


-kh K sinh ky - k cosh ky 


aD 
tt 
h = 
Pom Xo Fs ) £ L(k)e K cosh kh - k sinh kh 


oO 


cos kx dk 


for |x | >@, y =O as required in Section 3.32 may be determined 


by considering the complex contour integral 


K sinh ky - k coshky _-kh ikx ee ikx 
eee) eGeh Khe k siahkh © - di = 6 F(k)e ons 


where the variable k is considered complex, temporarily. The 
integrand of 255 is an even function of x so that we may use 
| x| in the following. It is understood that h >0O. The contour 
of integration consists of R >Re k>O, #£=indented above the 


Singularity defined by the positive real value of K Satisfying 


K coshkK h =- K sinhK h= O 
O fe) fe) 


for given K = hae and h>O, the part of the imaginary axis 


R>iImk>O, and the arc |k| =R, Re k >O, Imk>oO. Here 
Rn 7 Ki this is no restriction since we finally shall let Ro. 


integration around this contour yields 


e e Mi Lkx 
sr(eyer™ ane f Fleet ak + f Fleet ™ ak + f Flee 7 *ak=0, 
—_ —P n 


~ 


with symbols indicating segments of the path. With the contour as 


aefaned, 
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imGc\asedtess. 6 stems 


an 


Therefore after letting R-—->q@ and the radius of the indentation > O, 
we have 


ve MGcyes dia niR + /f F(k)e rk !* | dk , 
i o } 


and, after taking the real part, 


aD 
O56 %9¥ 9h) = Re miR. + ff [Re ir(ik)]e7*!*| dk 
O 


where k is now a real variable and R. is the residue evaluated. 
below. 


The residue R is 


R. = Lim (c-K, Le" K_ sinh ky - k coshky  ik|x| _ 


k=K, K cosh kh = k sinh kh 


The limiting operation is executed after the substitution 


= i h h 
K K sinh KS /cosh Ky 


and manipulation of the terms 


Go - K,) 


(k-K.) cosh Kh 
k cosh kh = k sinh kh 


K sinh (k-K )h+ ( * 
Oo O 


k-K.) Sinh kh cosh K oh 


With this substitution, the residue simplifies to 


“eR ik, |x| 
Ro = L(K )e E K cosh K. (h-y) : 


(k-K.) 


j = keK sinh Kh cosh K h 
K=K, K sinh (k K,)n+ ( =) _ 





Lh 


Ovegusencer evaluating the Gdamic. 


-K h 
7 O en : 
Re niR, = K LK, de 2K h + sinh 2K fh cosh K | (h-y) sin KS | x | . 
We find 
“Kh K em em K es .es 


sre O 


N 
i Oe Di a ene, eae 
eo ~ (2m-1)! cosh Kh bi 26) (2s-1)! cosh Koh 


by using the definition of L(k) in Section 3.31(c) and substituting 


K 


Ko Sinh K h/cosh Koh as well as evaluating the expression for 


k=K_. 


0 


With this expression we may write 


cosh K (h-y) 


Re miR. = - EB, (Kja,K hb) cosh Kh sink |x| 
where 
K em em N = K es cs 
Se eT O r nye 2n+l 0 . 
cemeast St) = 2K h + sinh2K nh | (2m-1)(* © (25-1)! 


The result for Re miR. is identified as the asymptotic behavior 


of Pon Xo ¥ oh) Lor | x | > om. On the surface y= 0, 


P56 %10yh) ~ -§E sin K | x | as |x| om . 


om 


We note that Eon >0O as hwryo@. 
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A-2 Asymptotic Behavior of ~_(K,h). 


The asymptotic behavior of the function 


'€ B) 
: cosh k(h-y ) cos kx 
Da) ENE? Sree ara = eee ae 


for | x | > om as required in Section 3.34 is found by following the 
procedure of Appendix A-l. 
The contour integral to be considered is 


$ cosh k( hey) otk |x| 


Rens ainemiccetnineie 2 


The contour is taken to be the same as in Appendix A-l. After fol- 


lowing the procedure outlined there, one finds 


(8 8) 
9 .(K,h) = Re niR. ef [Re iF Galea eK lI dk . 
.e) 


Thus the Re miR - gives the first term in the asymptotic expression 


for p (K,h) as |x| > om; that is, 


- ant 
9 ,(K,h) 2K h + sinh 2K h 


cosh K hcoshK (h-y) sin K, |x|]. 
On the surface y=0, as |x| >o@, 


ow ( } 
9,(Kyh) 6 E(K2) sin Ky |x| 


where 


2n (cosh K n)® 
EK h) = 2K h + sinh 2K h : 


We note that E (Koh) ~nt as h?®o@. 
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A=-3 Expanded Forms Useful in Numerical Calculation of 


Certain Definite Integrals. 


A-3.1 The Coeciircients of P56 6) s Qo,66) in Expansion 
Equations for the Method of Integration, 


Coefficients a and b . are defined for use 
em,k 


2m, j 


in solving the expansion equations by the method of integration 


(Section 5.12). 

















n/2 
4 f,(9) sinko de = 2 om. k 
(2m+k-1) 
2 
k odd: ene! k odd: k = n= Ll: O 
um = Ke rs 
k=1A4¢m: (-1) r/4 
k = 2m-1 41: + 1/4 
5 ; 
- Stone otherwise: 0 
(2m+k-2) 
k  2m(2m-2)(-1) 
k even: k = 2m: - 1/4 k even: —5—— ( 2 ( 5 
k -1l (2m-1) - k 
k # 2m: O ( 
n/2 
J £69) cos je dao = Dom, 5 | 
Deere = 
em 5 m e 1-) 
>) iis bm? - 42 + oy (-l) 
J Ani ere la 
5 odd: 1+ 
j+em-1 a 
5 even: as 
fe) 
7 ey a 
j+em-1 if 
2 odd em=-1+j 
Item even: @) 
. 7 
gee a AE)" a{* 4 
ie Coime yee 1-3 oe re 
| 


(zero is treated as an even number.) 
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R= 522 Two Special Forms. 


The following definitions are used for convenience 


of notation in later integral expansions. 


t/2 
I,(k,n) — yf sinkO sinnoe de 
O 
O if k or n are zero. 
Es 
=<+ 7 forn=k = 10 
fOr ovher ik, n: 
(n+k) even: O 
(n+k) odd: 
if k even: (-1)6atk+3)/2 Cae aan) 
a 2 
if ke odd: (a2) SBPEHLIZ2A nn? k®) . 
n/e 
I, (n,k) =f cos n@ sink@ do 
O 


O for k= 0 


for other K,n: 


If (n+k) even, and if (n+k)/2 even, then 
if k even: O 
Kk odd: 1/(k-n) 
if (n+k)/2 odd, then 
if teteven:. Cee) 


k odd : 1/(k+n) 





wf (n+k) odd: 
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A-4.5 Expressions in the Right-Hand Side of the Gey), 
E Pom 


Expansion Equations. 


t/2 sin ke fe) P I(k,p) 
js Yon 6 O99) = TT >: (-7)Pt1 a = ’ 
fe) cos je p=0 Pp. I, (jsp) 
/2 sin ke Tack 

ff  ¥,,68s"/2) d@ = n sin6 : 
fe) cos jo T,(js1) 


A-3.4 Alternate Forms for P6558), Yea 65s0)- 


The form of the quantities Yoq 6598) and ~ (648) 
used in the Qo,,65) expansion equation (Section 5.12) and in the 
calculation of M(6,@) (Section 5.22), respectively, are derived 


from the Cauchy Principal Value integral 


wo ky +ikx 
fat _ _-K(y-ix) 
d = ay ls EB, kK(y-ix) J 
-K aD phe tne 
=> * (cos Kx+ i sinkx] [y + log. r - i0 + Reel 
_ n=. : 


Here E, (-iz) designates the exponential integral of a complex 


argument [Tables of the Exponential Integral for Complex Arguments, 
National Bureau of Standards Applied Mathematics Series 51, May 1958]. 
Euler’s constant 0.57721...is designated Y and ©'= arctan x/¥ 

The real part of this expression represents a source in the free sur- 
face of fluid of infinite depth (Section 4.33) and the imaginary part 
represents the conjugate stream function. Contour integration leads 
to the real and imaginary parts given in Section 4.2, but the right- 


hand side of the above expression is one suitable to numerical 
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calculations. Evaluating this, 


6518) = aes oe [Q(0) cos (6sin®O) + S(@) sin(6 sind) ] 


¥.,(80) = Br oo ea aconi sini (seine ne Coco ca ntone 





where 
ae 
S(0) = 9+ 2 —=— sinno 
nin 
n=l 
oO ,n 
Q(0) = y + log 6 + ny aTn COS NO - 


It is noted that 


S(n/f2) = n/2 = Si (5) 


Q(n/2) = Ci (6), 


where Si(6) and Ci(&) are the sine and cosine integrals. 


A-3.5 Expressions in the Right-Hand Side of the do 6 6) 


Expansion Equation. 


The following integrals appear in Section 5.le. 


n/2 t/e rey sin ko 
ip ¥gq65:9) sine ao =f > Q(0)sin (6 sin9)sino d@ - 
0 O cos jo 
n/e Sin ko 
~-f 67° °°S® 5(9)cos (6 sin @)sing ado 
O cos je 
> ptl 6 T2E1P) ar 2 p+l 6 oe 
=(y + log 5) ?. (-1) ae . + 7 Alle eG) ir . ‘ 
p=0 I,(j,p) n=l p=0 I, (j,p+n) 
T,(kyp) 


(6 8) 
+ 7 ene 2. 
p=0 Po I,(Jjop) 
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The quantities not previously defined in this appendix are 


T(jyp) = 5 CAC|p-d]) + AC{ped|)] 





where 
2 
L ;t 
= (5) for $ 5c] R= Our 
A(|x|) = + [-1 4 (-1)*/2) for |x| even, 
x 
1. 1 (x+1)/2 
so er ei oe, for |x| odd; 
x 
and 
I, (k,yp) = [B(|p+k]) + (signum k-p) B( |p-k|)J, k #0 
O for k=0 
where 
O for |x|=0O. 
tt (x+2)/2 
B( |x|) = sy (-2) for |x| even, 
A (21) 72 for |x| odd. 
x 
n/2 sin ke n/2 sin ke 
f  y¥..(6,2/2)sino dg = ncos6 f sing de 
O . cos joe Oo cos je 
© ,.-x n/2 sin k@ 
+fos f sine ax de 
o0 6 +X fe) cos je 
I,(k,1) 


= - [(n/2 + Sid)cos& - Ci 6 sin 6] 


I, (452) 





sone 


This form is found directly from the expanded form of ¥o,66s0/2) 
given in Section 5.11, and the identity 

© ..-x 

vs > «dx = [n/f2 - Si (6)]cos(6) + Ci (6) sin (6) 

o 6 +x 
(Lamb, Hydrodynamics, Art 244(37)]1 Here Si(6) and Ci(6) are 
the sine and cosine integrals and are available in tabulated or 
series form. Alternately, the same result for Yoq65st/2) may be 
found in Appendix A-3.4 above where a series form suitable for com- 


putation is given. 
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APPENDIX B 


The Finite Deflection of Circular Plates with Initial Tension. 


B-4 


B-5 


B-6 


B-/ 


B-8 


Pitroduction, 
Governing Differential Equations. 
Series Assumed for Solution. 
Evaluation of Coefficients don in Terms of b 
and Definition of the Function Cy 
Evaluation of Coefficients Bon in Terms of b 
and Definition of the Function Coe 
Evaluation of b, by the Boundary Condition. 
Graphical Solution. 
Special Cases. 

a. Zero Initial Strain. (The Solution of Hencky). 


b. Lightly Loaded Plates.’ 


¢. Comparison with the Membrane Solution. 





Boo), 


APPENDIX B 


The Finite Deflection of Circular Plates with Initial Tension. 


Be] MNeroaguction. 


Consider a circular plate of homogeneous material obeying 
Hooke’s law which is initially flat until uniformly pressure-loaded 
on one side. The edge is clamped. The only widely-published solu- 
tion of the simplified equations of equilibrium including a given 
initial-tension boundary condition is the membrane solution. The 
membrane solution assumes vanishing bending stiffness, small de- 
flections, and that the initial tension is large compared to de- 


veloped material stresses. 


In order to maximize pressure sensitivity of a diaphragm- 
type pressure gauge, one is led to thin diaphragms of vanishing 
bending stiffness and to reduce the initial tension to the minimun. 
Application of the membrane solution is not wholly satisfactory 
when initial tension is small because material stresses become ap- 
preciable at moderate loads, and because deflections may not be in- 
finitesimal and because the solution gives no hint of departure from 
linearity. The exact solution of Hencky assumes zero initial tension 
and the infinite deflection derivative at zero load is not physically 


appropriate to stable diaphragms of finite pressure sensitivity. 


This appendix presents an exact solution to the equations of 
equilibrium for this type plate subject to arbitrary initial tension. 
For zero tension the solution of Hencky is recovered. For large ini- 


tial tension the membrane solution is approached asymptotically. For 
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lightly loaded plates the solution always tends to that of the mem- 
brane. The initial deflection derivative is limited by the initial 
tension, physically appropriate to the case of a lightly loaded, 


thin, stable diaphragm. 


The radial and tangential stresses, as well as deflection, 


are given in this solution. 


The membrane solution (large initial tension) is valid only 
for infinitesimal deflections. This solution is not so limited. 
The deflection pressure relation of the membrane solution is linear. 
This solution reveals the nonlinear solution, tending toward a 
Hencky-type solution as the developed stresses increase relative to 


initial stresses. 


Bee Governing Differential Equations. 


The equations of equilibrium and stress compatibility for an 
initially flat circular plate with clamped edges were given by 
t 
A. Foppl [1907] and are repeated by H. Hencky [1915] and S. Way [1934] 


among others. 


We repeat the equations expressing horizontal equilibrium, 


da as 
dr (or) 2 Shes Os 


and compatibility, 


2 
d E ,dad _ 
Ties Oat J, ) + 5 C=) = O'%4 


using oe and o, as radial and tangential stresses and 6 as de- 


flection. 
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The equation expressing vertical equilibriun, 


5 
bP t da (lh a ¢, o6)) _ pr ede) 
1 Bip ae ae ne ae ee 
me a 
is simplified in the case t?<< Gs 
"ee ee a6 


2a ruar ~ 
This is the case of vanishing bending stiffness. For extended argu- 


tt tt 
ment see, for example, Foppl or Hencky [loc cit] or A. Foppl [1924]. 


We deal here only with the case of vanishing bending stiffness. 


Be 3 Series Assumed for Solution. 


A power series form is assumed for the dimensionless radial 
stress ratio (o/E). The corresponding series for the tangential 
stress ratio (6, /E) is derived with aid of the equation expressing 
horizontal equilibrium. That differential equation is thereby satis- 


fied. 


The stress is made non-dimensional by Young’s Modulus E; the 
radius r, by the plate radius a. The dimensionless coefficients 
bon are to be determined and Cy is a dimensionless function of 


the uniform pressure load found later to be 


2 
3 _ 1 ;pa 
c,° = gy CEP 


From horizontal equilibrium, the series form for 5, /E is 
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tp, a(n), & 
Be age oe ie 
5 - on 
= C Con-ae 1) <b (— ° 
40 =O en-a 


A series form is assumed for the dimensionless deflection 


(6/a) compatible with the condition 65 =O at r= a. 


bees [Pan - Fag Om 
sae a - a (=) ‘ 
a 2 el en nz0 Cn ‘a 
The dimensionless coefficients a are to be determined 


en 


and c is a dimensionless function of the uniform pressure load, 


2 
later found to be 


The center deflection at r=0 is 


aD 
6(0) = (ac,) > as, ° 
n=0 
B-4 Evaluation of Coefficients don in Terms of b, and 


Definition of the Function Cy 


With aid of the equations of compatibility and vertical equi- 
librium we are able to express Bo, in terms of eS after making 


a suitable choice of definition of Cy) 


The nonlinear equation of compatibility in the form 


2 P 
(ory ce aay, ty ‘ 1 (or ab, =O 
E Y dr ‘E E 2 ‘rE ar 


with the equation of vertical equilibrium 





EO7, 


Q 


95 _ = dspr 
Ear 2 Et 
takes the form 
a a Je Oe af pr,* 
=m as ne mo See) Oe 


The first term of this expression is rewritten as the product of 


two appropriate series, using for one factor 


q fe) e en fe) en 
r a = =c, © | 2 eaD + Ae bo, (2) | 


r.* fe 9) r.cnee 
cy cae OO iain b5 6S) ’ 


so that 


ane 2r- a - en c = en-2 1 > 
Cy =) Pa | by: (2n+2)(2n)b, (=) | + z Ce) =O. 


= nat 


It 1s appropriate to define 


4 
2) ee a 
eC e 


Cc 


in order to leave a non-dimensional expression which is used below 
to generate recursion formulae for Don? 


en=2 


2 
@ en a 
| 2 Pane | [z (2n+2)(2n) v5, (E) |+ ee. 


n=1 


This series is expanded and coefficients of powers of r/a are set 


to zero individually, the first reads conveniently 


8 bob. + 8=0 or bob. e 1 =O os 
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The succeeding equations for the next higher powers of r/a are 


a 2 
Dd) d + ebob. : 


Se roe etc. 
Oo 


2 
bbb. a 8b),b A a= 


2 


Solving these in succession, all succeeding coefficients bo, are 


expressible in terms of La The first five are 


b , to be determined from boundary condition 


O 
b, = - 1/b,* 
b, = - 2/3b. 
be = = 13/18» © 
be = - 17/18b,7* 
b107 ~ 37/27," ° 
Be5D Evaluation of Coefficients a5, in Terms of dD, and 


Definition of the Function C5 ° 


With use of the equation expressing vertical equilibrium we 
are able to find all coefficients Bon in terms of bo after making 


a suitable definition of the function C5 ; 


The form of the equation of equilibrium 


Q 


= d5 , i pr _ 
Reece: to cane 


with the series for o/E after using 


AG eo) pent ae ron 
aa = = 2 Oh iar es =o ccs a eee as (=) ) 
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is 
aD . 2n @® 7 en 1 pr 
ro PD 2 © PE Coed) D> (nel) aj, (+ +c =O. 
Ieee n=0 en'a Lay, en’a 2 Et 
Substitute the previously selected definition of Cy and define 
o =i Pal_ 3/pa 
“2 3-8 ttc, > APE 


8 


in order to leave a dimensionless expression which can be solved 


Successively for a 


Dae 
@ r. @ ‘ on 
Pinder) Ga), Kaede eG Rae sO ue 
n=O n=0 


When this series is expanded and the coefficients of powers of r/a 


set equal to zero individually, one has a series of equations 


0 


abel 
oO 


O= 2a5b + a ios 


O= 3a},d. + 2a5d., + a Py, ~ ‘GCs 


Suesertuting for oon in terms of dD, and solving successive 


equations yield for the first six Aon 


ape rat 1/b, 
h 
a, = + 1/2b. 
7, 
2 LO 
ag = + 55/720 | 
ag = + 105/90b ~~” 


1.6 
ake +- 205/108b | ° 
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B-6 Evaluation of b > by the Boundary Condition. 


The unit strain e¢ in terms of the in-plane stresses is 


The edge r=a of the plate is rigidly clamped. If the plate 
is not initally strained, ¢€¢ =O at re=a and must remain zero be- 
cause of clamp action. If the plate is initally strained and then 
clamped at re=a, the strain ec = Eo at re=a and by clamp action 
must remain Eo? The boundary condition is 

S. - po = Be 4 ; 


a & 


initial strain. 


Ee 
@) 


The series form for 5, /E and o /E at re=a are 


e9) 
oJ = cjE 7a (2ént+1)b,, 
r=a n=O 
e9) 
oJ = cyt 2 Le ° 
r=a n=0 


The boundary condition in series form is 


cD 
aie = pbs, = Ee /ey 


or, expanded, 


oO 





anonk 


For given material and load this equation can be solved fora 


numerical value of by and the problem is solved. 


B-7 Graphical Solution. 


A convenient procedure for evaluating dG from the boundary 


condition equation and the deflection coefficient K = 2 is a 


a 
n=0 © 


graphical one, 


The left-hand side of the boundary condition equation, includ- 
ing terms up to the n-power in b , is designated f fe, so that 


with sufficiently large n and db, ? l it is nearly true that 


f 6b.) = e fe (1 - p) . 


Values of bo were assumed in the range 1< be < 10 and n was 
taken as 17. The numerical values of f,5(b,) are plotted in the 
figure to enable graphical solution of the inverse problem: given an 
initial strain E the quantity Efe, (1 - p) is computed and the 


corresponding dO may be determined graphically. 


The center deflection is 
= pa 
6(0) = (ac5) 2 Fan aN ate 


where we have introduced K, defined by 





\@ 8) 
Ih i gz, 22 (de) 205 
Ka 9) a = + + ee, + ee © 
neofee 3G 2p 9p! 72b. 60b,.7 108b_* 


This summation has been computed for the assumed values of Boe 1<b <0, 


using these six terms, and the numerical value of K is presented 
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THE FUNCTIONS 4. (b,) AND K(b) 


FIGURE 4Z 





163 


graphically. As a matter of convenience, the two figures 


£i5(b,) vs do 


are graphically combined to eliminate bo and present 
K vs é (fc, (1-p) . 


With this graph, one enters with known load and initial strain to 


find a value for K and with this the center deflection is 


B-8 Special Cases. 


ae Zero Initial Strain. (The Solution of Hencky). 

If the initial strain is zero, this solution be- 
comes that of Hencky. The series expression f(b.) is equal to zero. 
The corresponding value of K is evaluated for the root dD, of 
f(b.) =O. Our computed value using £55(b,) differs by less than 
one percent from Hencky’s value K = 0.662 and the difference is no 


doubt entirely due to numerical procedures. 
Die Lightly Loaded Plates. 


This special case is appropriate for considering 
the initial pressure sensitivity-and for comparison with the zero- 
initial tension case which predicts an infinite initial-deflection 


derivative suggestive of an unstable initial position. 
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As the pressure load approaches zero, Cc, ap- 
proaches zero, and f(b), proportional to 1/c,, is increasingly 
large. The function f(b.) approaches its asymptote, d, 


3- 3 
f(b, )) = b-oe > - oo. yb . 


o l- fe) 
be large ii Le 


For large values of Dos the value of K is 


We 


cD 
_ 0 iL 1 
K = > So) b 7 a pe ae eee bd 
n= fe) 2b. fe) 


ne 


One may substitute for f(b) its equivalent in 
terms of the initial strain, E /c,(1-p)s but it is more revealing 


to make another substitution as follows. 


Under the assumed conditions, bo is large and 
therefore all Don for n>O are negligible. Therefore, at the 


edge r= &, oe = o, and the edge boundary condition may be written 


o, - Bo, = o(l-n) = wo r 
r=a 


The initial tension T = ot, so that one may use for f(b.) an e@X- 
pression in terms of T: f(o)) = E /c (1p) = T/c,tE . With this the 
center deflection coefficient K is c, tE/T and the center deflec- 
tion reduces to 


pa* 
SCO) f= Lt 


which is exactly the membrane solution. 
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Thus the initial deflection of a lightly loaded plate is con- 
trolled by the initial tension. It is not considered appropriate to 
assume zero initial tension for a stable diaphragm of vanishing bend- 
ing stiffness; the initial deflection derivative is bounded and equiv- 


alent to the membrane result am /kn , 


Co Comparison with the Membrane Solution. 
The membrane solution may be recovered here by 
observing that a particular solution to the equation expressing hori- 


zontal equilibrium is 


o. =o = Oo = constant =T/t . 


The equation of compatibility will be satisfied only for in- 
finitesimal deflection in this case, so that Giey aa is negligible. 
The vertical equilibrium equation is equivalent to the membrane equa- 
tion, 


96 — ee) yoy ers 


In terms of our series solution, Bon are zero for n>?O and the 


result is given in section b above. 


The complete series solution given here is not limited to in- 
finitesimal deflection and shows that as the pressure load (deflec- 
tion) increases, C4 increases, f(b.) decreases and the deflection 
tends to that of Hencky [f(b ) = 0] as the developed stresses be- 


come increasingly large compared with the initial stress. 
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APPENDIX C 


The Attenuation of Oscillatory Pressure Fluctuations Travelling 


in a Viscous Fluid Through a Connecting Tube to an Instrument 


Cavity with Resilient Walls. 


C-4 


C-5 


C-6 
EY, 
C-8 
C-9 


C-10 


Description. 
The Navier-Stokes (Momemtum) Equations. 
The Equation of Continuity. 


The Equation of State and Effect on Tube 
Wall Elasticity. 


Boundary Conditions and the Effect of 
Diaphragm Elasticity. 


Solution of the Governing Equation. 

The Short Tube Approximations. 

Approximation Valid for Small Values of Wey) 
Approximation Valid for Moderate Damping Values. 


Approximation Valid for Small Damping Values. 
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APPENDIX C 


The Attenuation of Oscillatory Pressure Fluctuations Travelling 
in a Viscous Fluid Through a Connecting Tube to an Instrument 
Cavity with Resilient Walls. 

C-1 Description. 

The required diameter of a diaphragm-type pressure gauge may 
be controlled by the sensitivity desired. This diameter may exceed 
the allowable diameter on the surface where pressures are to be mea- 
sured. The apparent diameter can be reduced by communicating the 
pressure through a port of allowable diameter to a cavity whose one 


wall is the required diaphragm. 


The apparent diaphragm diameter is reduced, but there may be 
a penalty in degradation of the dynamic response of the system tube- 


cavity-diaphragm compared to that of the diaphragm alone. 


This appendix reports an analysis of the attenuation of oscile- 
latory pressure fluctuations travelling ina viscous fluid from an 
entry port through a connecting tube to an instrument cavity. The 
instrument cavity is expected to have a resilient wall, the diaphragm 
of a pressure gauge in our application. Related work is known in the 


literature but is usually concerned with long tube lengths, compres- 


sible gases, or does not include the effect of a resilient wall. 


The present analysis concentrates on the arrangements expected 
in measurements of small pressure fluctuations on the surface of 
small bodies oscillating in a fluid. The slight compressibility of 


a fluid and the high flexibility of a sensitive diaphragm contrast 
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with the more usual situation of gas compressibility and relative 
cavity rigidity reported in acoustics. Diaphragm elasticity plays 
an important role. It may in fact be the limiting factor in dynamic 


measurements of this type. 


C-2 The Momentum (Navier-Stokes) Equation Yields Particle 
Displacement and Volume Rate of Flow. 
The momentum equation for a fluid element in the tube is the 
Navier-Stokes equation for the axial velocity component w. We as- 


sume radial symmetry. 
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2 2 
Ow Ow _ g Op Qow 1 aw Qow 
ge az yea Ot ee a ee 
or Oz 


Assume that dw/éz is negligible compared to d@w/dt and aw/ar. 


Assume also that the particle displacement is harmonic in time. Then 


6(r,t) ocr) ; 
w(r)cIM® o> = jw6(r)e 2% : 
2 


SY 


ees: 1 (2 We ae + o8) _ & 9p 





9 
t aré r or p Oz 
2 
d“6 wo Onno) 


2 radar 


A particular solution of this equation, assuming p not a function 


Of er, is 


= (SP) o/y- 
6 (55 8/H Pp : 
A homogeneous solution of this equation, finite at r=0, is 


Oe J Cxr) 


Race ju/Vv . 


Combining these solutions with a constant c, 
6 = (dp/dz)g/wip+ co J. (kr) 
a ae i O ; 


The constant is evaluated by the no-slip condition that 6 =O for 


r= 4. 





— 
— 


(22) ¢ 3. (kr) 
a al eal = ® 
2 | re | 
sa 
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It is noted that the radial derivative of particle displace- 


ment, 06/ér, vanishes at the centerline r=0 as it should. 


The volume rate of flow through a given section is found by 
integrating the area-velocity product over the tube radius at that 


section. 


. 96 jong (ap/dz) Jo (kr) cdr 
qaV = s> enr ce Jr ax - TE) | 
a 
vefs av 
O 
2 ka 
_ J2ng(dp/ez) a +-— } kr Jo(er) (kr) 
fe kJ (ka) O 


jna“g(ap/dz) E ; ed, (ka) 


pw kad (ka) 


For simplicity we use notation 


le 23, (ka)/kad (ka) 
V = - BCdp/dz) 

Vd 
B = na g(1-3)/j pw ° 


It may be observed that for small (ka) 
2 
1-7 > (ka) /8 ° 


Substituting this in the expression for the volume rate of flow, 


V=e- na ‘g(8p/02)/8 p 4 


the known volume-rate of Poiseuille pipe flow in which inertia forces 
are negligible. That is, for w/v vanishingly small this solution re- 


covers the Poiseuille solution. 
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C-3 The Equation of Continuity. 


The equation of continuity of fluid in a thin section of a tube 


of area A and flowing with radial symmetry is 


“fe otgea ai) $ 


We will substitute the previous result 
V = wA = = B(dp/dz) 


and change variables from the actual local pressure p to the frac- 
tional overpressure 9». The overpressure 4, or pressure signal, is 
normalized with respect to the amplitude ise of the sinusoidal pres- 
sure assumed to exist at the tube entry. The forcing pressure Py 


is assumed to be independent of r. 


9 = (p = p)/P, 
Op _ Oo” 
az ~ Po az 


aph 2 a95 _ 
Se + a5 [ - php, gg) 20 - 
c-4 The Equation of State and Effect of Tube Wall Elasticity. 


Introduce now the equation of state of the fluid in the tube 
in the form of a relation Pp P) . For the small pressure fluctuations 
expected in this work the density-pressure relation may be expressed 


by means of a constant bulk modulus of compression B of the liquid: 


aa 
Pp 
P= Poo * a? 


Analysis from this point concerns nearly incompressible liquids. For 
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gases, a similar path would be followed but different equations of 


state are necessary (A. Iberall [1950]). 


The elasticity of the tube walls may not be negligible in com- 
parison with the slight compressibility of the liquid. In this case 
the cross-section area A is a function of pressure and might be 
written as approximately 


2 7 *o) 


A = AG eat 7 


where E = Young’s Modulus of wall material, 


f = a geometric function of wall thickness and nominal 


form. For a thin-wall circular tube f = a/t, 


A, = cross-section area at p= Po: 


[For a table of f for ordinary shapes see I.P. Ginsberg, Prikladnaya 


Gidrogazodinamika, Leningrad University, 1958, p. 106.] 


Combining these results to define an apparent bulk modulus in the 


tube, we find 


P- Py 16 ee ee 
ee ee 
- P 

~ PoAy St j B =) : 


iW aaa 3 
B B, E 


Whether this apparent modulus B differs from BL is a measure of 


the effect of tube wall elasticity. Using the above expression for 
(pA) » we obtain 


ot 2B 


a] & 
e 
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With this result, the equation of continuity is 





Poo a ac 
B ae ae ee 


Here we have required that the product PP be independent of z. 
It is noted this product is independent of density, so that we have 
not assumed density independent of 2z but we have assumed that the 


kinematic viscosity is sensibly independent of z. 


C=5 Boundary Conditions and the Effect of Diaphragm Elasticity. 
The boundary conditions are these. 


(i) The pressure fluctuation at the tube entry is specified. 


gOS 


p= 
= Zo 


(ii) The mass rate of discharge from the tube equals the 
accession of liquid in the cavity. This augment to the 
contained fluid is accommodated by elasticity of the dia- 
phragm and compressibility. 


The second boundary condition is expressed analytically as follows. 


The fluid mass in the cavity is p* so that the mass aug- 


ment is av 
do 


O 
[ye oS | 
From the previous relation ple) we derive 


0 


1 


ate 
"Oo 
a 


To determine the rate of change of instrument cavity volume we in- 


troduce the volume sensitivity, s, of the diaphragn, 
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The volume sensitivity is geometrically related to the pressure sen- 
Sitivity o of the diaphragm, a measure of the change in center de- 
flection per unit change in pressure. For a membrane diaphragm de- 
flecting to a paraboloid of revolution, 

Site EVE. : 


With the volume sensitivity defined as above, we find 


The Mass influx into the cavity is 





The expression (Vo + SB, ) represents an equivalent cavity volume 
V corrected for the resilient diaphragm. The term Poof 31 is 
equal to a/oX, the squared velocity of acoustic propagation in 


the volume, With these notations the mass influx 


oD T/c* 
a1 = Jw piev/c 


Pre 


is equated to the mass flux at the end of the tube, 


0 OY 
va = - p8 <P} =- ppp. =] ‘. 
Vs a=L re Oz wee a Oo 02 ay 


The second boundary condition is 


ez 2 : 
Z=1 pBe 





a2 


where 


b = jwg T/oBe“ ; 


C-6 The Solution of the Governing Equation. 

The differential equation from Section 3, after combining the 
equations expressing momentum conservation, continuity, and the ap- 
parent equation of state, is 

2 
0p _ AOR ap 
a2° Ree ot 
The solution is 
AZ, “OZ y _ dwt 


where 2 
oa = juA. p,/pBB : 


The desired form of the result is 0/®, which is 


C C 
[= <% me ey! i [= + Dindiiaes 
2 2 


thes ratio c,/C5 is evaluated by means of the boundary condition of 


continuity at the cavity, 8 + bo = 0 at 2ae=tL, so that 


oe (ashe 
com “ae 


Finally, evaluating the pressure-Signal ratio at z~a~=lL, 


ce (aL) 
® ~ (aL) cosh (aL) + (bL) sinh (al 


The ratio 0, /?, is the pressure fluctuation in the cavity per unit 
fluctuation at the tube entry. This ratio is called P;/P, in the 


text. The products aL and bL are as follows. 
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Pe eck Cee ree ee 
“, Vlie- J 


Pe 
Here c, , the apparent sound velocity in the tube is Po/ Bes 
determined by the apparent modulus B rather than the fluid mod- 
ulus alone as in the cavity. The ratio w/c, is the fluctuation 


wavelength in the tube, i.e. 





L 1 
aL=je >; ——. 
AVI-4g 
Similarly, for bL we have 
, Tat Beet 
bL = - w > Too ° 
cA 


OC. yan eee 
The quantity c A/VL is the square of a characteristic frequency, 





we = o“A /TL : 
n 
so that 
2 
- Ww 1 
biv= : > Tsa° 
n 


It should be noted that aL and bL are, in general, not 
real variables, as J is the ratio of a complex quantity and Bessel 


functions of a complex argument. 

The natural frequency we is lowered by increasing diaphragm 
sensitivity through the effective volume V. 
C~7 The Short Tube Approximation. 


If the tube length L is a small part of the pressure fluc- 


tuation wavelength in the tube, the product aL is small and 





a aL 
ie aL) cosh (aL) + (bL) sinh (aL 
~ al 1 
~ eb i 
ties cs ede nen 
W RE J 


This is the form discussed in the text. 


C-8 Approximation Valid for Small Values of Wen 
The value of 


; Pees 
~ ka J (ka) 


for small values of the argument 





(eae a jwaS/v 
1s approximately 
2 
J=Zzle (ka) 
s0 that 
1 ~ j8v 
i et ia a ee: 
wa. 


The response ratio. for this case is 


P. 
1 ab Cue 
— — a —————————————— = a 


This result is identical with the one obtained by assuming 
uniform Poiseuille flow and constant pressure gradient. If the 


Poiseuille flow mass rate, 


ed 
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h 
_ ~ ma 0 
Qapv=- yh (SY) , 


with constant pressure gradient, 


the governing ordinary differential equation is 


dp. ; 
at * Py 





v = p(t) A 


Assuming periodic flow yields the result given above. 


The time constant t+ is increased by increasing diaphragm 


sensitivity or decreasing the entry port diameter, for 


c= BILV/na to” , 


CH9 Approximation Valid for Moderate Damping Values. 


ine identity 


1 J (ka) 


ia ~ Jo (ka) 





is, for moderate values of the argument, approximately 


, 8 
(wa*/J) 


Wy 


The corresponding response ratio is 


7 ak 
= 9 
O ea) + ju rages 


‘oto 
He 





eh 


where we have used the notation 


This form is motivated by analogy to that for the response of a 
moderately damped spring=-mass system. This system responds with 
90 degree phase lag at a dimensionless frequency u=1.0 with an 


amplitude ratio 1/2¢. 


This form is also analogous to the classical Helmholtz 
Resonator results, but the present more-detailed analysis predicts 
a lower resonant frequency of 90 degree phase lag and generally 


greater phase lag at frequencies near resonance. 


To clarify the differences in the results, compare the above 
with the following. We recall from section 8 the ordinary differ- 


ential equation 


dp. 


SR ee oe 


Multiplying by the tube area makes more evident that this is a 
force balance of pressure force and Poiseuille tube-wall friction 


considered to be constant along the tube length: 


et ee ete © 7 
Me re a ena-L 3 = INC Pp, ) ; 
mag 


where (4vQ/na’g) is recognized as the constant Poiseuille-flow 
tube=wall friction per unit area. Add to this force balance an 


inertia term 
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2 
ap. 

oe ena : mee A MA a c“A , 
feo FP ee oe byes vee 


The steady-state solution, assuming periodic pressure fluc- 


tuations of frequency w, is 








Pil 
a ines 
Cars + j22 
n n 
where 
w E o“n 
LV 


as before and 


ae w t/2 = by/a~w - 


This is the result given by Rayleigh and generally observed 
to predict resonant frequencies that are higher than measured values 


(M. Colby [1939] p. 220; J. Patterson [1952] p. 9). 


The more detailed analysis, admitting longitudinal variation 
of pressure gradient and a more general radial velocity profile, 
predicts acoustic resonance at a frequency about 14 percent lower. 
This further lowers the expected range of frequencies free of acous- 
tic effects. It is observed that in either case the effective in- 


strument cavity volume is used in calculating w,° 
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C-10 Approximation Valid for Cases with Small Damping. 
The identity 


1 J (ka) 
a si ~ J, Cka) 





tends to unity for large values of the argument. The response 


ratio is therefore 


> er 
Po ‘; ae 
~ 2 
WwW 
n 


The value of this approximation lies in the fact that it emphasizes 


the tendency toward undamped acoustic resonance. 
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